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ABSTRACT 

Cosmological hydrodynamic simulations robustly predict that high-redshift galaxy 
star formation histories (SFHs) are smoothly-rising and vary with mass only by a 
scale factor. We use our latest simulations to test whether this scenario can account 
for recent observations at z > 6 from WFC3/IR, NICMOS, and IRAC. Our simula- 
tions broadly reproduce the observed ultraviolet (UV) luminosity functions and stellar 
mass densities and their evolution at z = 6~8, all of which are nontrivial tests of the 
mean SFH. In agreement with observations, simulated galaxies possess blue UV con- 
tinua owing to young ages (50-150 Myr), low metallicities (O.I-O.SZq), and low dust 
columns {E{B — V)^ 0.05). Our predicted Balmer breaks at z = 7, while significant, 
are « 0.5 magnitudes weaker than observed even after accounting for nebular line 
emission, suggesting observational systematic errors and/or numerical resolution limi- 
tations. Observations imply a near-unity slope in the stellar mass-star formation rate 
relation at all z = 6-8, confirming the prediction that SFH shapes are invariant. Dust 
extinction suppresses the UV luminosity density by a factor of 2-3, with suppression 
increasing modestly to later times owing to increasing metallicities. Current surveys 
detect the majority of galaxies with stellar masses exceeding IO^Mq and few galax- 
ies less massive than 10^-^ Mq, implying that they probe no more than the brightest 
w 30% of the complete star formation and stellar mass densities at z ^ 6. Finally, 
we demonstrate that there is no confiict between smoothly-rising SFHs and recent 
clustering observations. This is because momentum-driven outflows suppress star for- 
mation in low-mass halos such that the fraction of halos hosting observable galaxies 
(the "occupancy") is 0.2-0.4 even though the star formation duty cycle is unity. This 
leads to many interesting predictions at z ^ 4, among them that (1) optically-selected 
and UV-selected samples largely overlap; (2) few galaxies exhibit significantly sup- 
pressed specific star formation rates; and (3) occupancy is constant or increasing with 
decreasing luminosity. These predictions are in tentative agreement with current ob- 
servations, but further analysis of existing and upcoming data sets is required in order 
to test them more thoroughly. 

Key words: cosmology: theory — galaxies: evolution — galaxies: formation — galax- 
ies: high-redshift — galaxies: photometry — galaxies: stellar content 



1 INTRODUCTION 

In the era of precision cosmology, the study of galaxy evo- 
lution is an initial value problem in which the goal is to 
characterize the processes that c onnect the initial co ndi- 
tions at 2 ~ 1090 (for example, iKomatsu et al.l [2010|) to 
galaxies, which had evidently begun forming by z = 10 (for 
example, IStark et al.ll2007l : iBouwens et~ai]|2009d : lYan et all 
120091 ) . The initial- value nature of the problem implies that 



detailed measurement of its earliest stages constitutes an 
immensely useful aid through its ability to anchor an under- 
standing of the later stages. Gas inflows, outflows, and other 
feedback processes impact high-redshift star formation his- 
tories (SFHs). Hence observational inferences regarding the 
SFHs at high redshift constrain star formation and feedback 
both at early times and at subsequent epochs. 

Cosmological hydrodynamic simulations robustly pre- 
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diet smoothly-rising SFHs JFinlator et al.1 120071 ') with a 
scale-invariant shape at z ^ 6. Therefore, any compari- 
son between simulation predictions and observations con- 
stitutes a (possibly indirect) test of this scenario. We have 
previously demonstrated considerable success in using cos- 
mological hydrodynamic simulati ons to interpr e t obs erva- 
tions of high-redshift galaxies. In lFinlator et al.l l|2006r ). we 
showed that simulations of a ACDM universe with star 
formation modulated by outflows with a constant wind 
speed and mass loading factor broadly reproduce the ob- 
served rest-frame UV luminosity function (LF) at z ~ 4. 
We confirmed previous numerical prediction s (for exam- 
ple, iDave et al.l I2OO0I : [Weinberg etall I2OO2I ) that stellar 
mass (M«) and star formation rate (SFR) correlate, find- 
ing a small scatter and a slope near unity. A quantitatively 
similar correlation has now been observed at all redshifts 
out to z = 7 IBrinchmann ct al. 2004; Elbaz et al. 20(J^; 
ai] [20074 iDaddi et all l2007l: ISalira et al.l I2067I: 



Noeske et 



Schiminovich et al. 120071: iPannella et al. l2009l : Labbe et al 



200E : iBothwell et al.ll2009l : lOliver et al.ll2010l : iMagdis et al 
2010I ). although observations at z < 2 tend to sug- 
gest a somewhat flatter slope of 0.7-0.9 that is asso- 
ciated with gala x y downsizing (IBrinchmann et al 



Elbaz et all l2007l : iNoeske et al 



12007a 



Daddi et al. 
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Sahrti et al] l2007l: ISchiniinovich et al.l l2007l : iBothwell et al.1 
2003 : iMagdis et aLlbOlOh . This correlation has two implica- 
tions: First, it implies that UV-selected galaxies are gener- 
ically th e most massive (star- forming) galaxies at a given 
redshift l|Nagamine et al.l 12004 : Finlator et al. 20061 ) unless 
the scatter is quite large (jWeinberg et al.l [20021 ). Second, a 
slope of unity supports a scenario in which SFHs have a 
scale-invariant shape such that, on average, galaxies' SFHs 
differ only by a scale factor. The SFR-M* relation there- 
fore constitutes a key prediction of galaxy evolution mod - 
els (for example, iBouche et al] l2009l : iDutton et alll2009t ). 
We additionally predicted that high-redshift galaxies should 
exhibit pronounced Balmer breaks because smooth SFHs 
naturally yield evolved stellar populations. With the ar- 
rival of rest-frame optical constraints from IRAC, this pre- 
dicti on has now received dramatic confirmation (for exam- 
ple, iBradlev et al.1 '2008': ' Chary. Stern fc Eiscnhar dj l2005l: 
Dow-Hygelund et al...2005;:lDunlop et al.ll2007.:.Egami et al. 



20051 : lEvles et al.ll2005l. l2007l: iLabbe et al.l I2OO6I: iLai et al 



2007^ : iLabbe et al' '2009': 'Stark ot al.' '2009': 'Verma et al 
20071 : IVan et al. 2005; Pentoricci ct al. 2009; Zheng et al. 
20091 ) ■ although the question of whether emission lines could 



be mi micking the observed Balmer bre aks remains unre- 
solved (jSchaerer fc de Barroj[2009l . l2010l ). 

In iDave et al.] (|2006l . hereafter DFO06), we tested a 
treatment for momentum-driven outflows from star-forming 
regions against the available constraints on the rest-frame 
UV and Lya LFs at z ~ 6, finding reasonable agree- 
me nt in both case s . This result was subsequently extended 
by iBouwens et al.l ([20071), who found that our simulation's 
prediction for the evolution of the UV LF normalization 
from z = 8 — )■ 4 was well-matched by observations. This is 
a nontrivial test of our predicted SFHs because constant or 
decaying SFHs would predict a different evolution of the UV 
LF than is observed. 

Finally, in iFinlator et al.l l|2007l ) , we showed that our 
simulations could account for the rest-frame UV-to-optical 
spectral energy distributions (SEDs) of 5 out of 6 ob- 



served galaxies between z = 5.5-6.5. We demonstrated that 
smoothly-rising SFHs are equally as plausible as more popu- 
lar constant or exponentially-decaying models. We then used 
these comparisons to explore the implied physical properties 
such as SFR, M«, and age, finding good agreement with pre- 
vious constraints although with significantly narrower pos- 
terior probability distributions owing to the tight priors im- 
posed by the physical correlations that our simulations pre- 
dict. 

While these studies are encouraging, their support for 
our SFH scenario is somewhat indirect because, until re- 
cently, observations at z ^ 6 did not clearly imply a unique 
SFH. For example, any galaxy evolution model that pre- 
dicted scale-invariant SFH shapes would predict an SFR- 
M« relation with a slope of unity. Similarly, any SFH shape 
can be tuned to reproduce the observed UV LF at one 
epoch, and bursty SFHs could be tuned to reproduce it at 
multiple epochs. The fundamental difficulty is that broad- 
band SEDs do not contain enough information to constrain 
the underlying SFHs of individual objects, even when an- 
alyzed throug h detailed SED-fitting stu dies (for example, 
IShaplev et aLi r2001 : Papovich et ahfeoOll ). 

Since the completion of these previous works, a num- 
ber of authors have used new data from the HST WFC3/IR 
and NICMOS cameras as well as Spitzer/IRAC to increase 
the number of objects at z ^ 6 with publishe d re s t-frame 
UV to optical con straints (iBouwens et al. 20 09bl. 2010a : 



Bunker et al.l 20091: ICastellano ot al. 2009; Finkelstein et al 



2009"; 'Gonzalez et al. 2 009i: iLabbe et a l. 2009; Labb e et al 



2010; McLure et 



al. 201S 



Oesch et al. 2010a, b; Ouchi et al 



2OO9I : lYan et al.ll2009l ). Using these new samples, we may 



now ask whether the mean SFH of reionization-epoch galax- 
ies is c onstrained e v en if individual SFHs are not. For ex- 
ample, IStark et al.l ||2009|) recently pointed out that the 
observed lack of evolution in the specific star formation 
rate (SSFR) of UV-selected galaxies from z =4-6 (see also 
lOverzier et al.ll2009l : lGonzalez et al.ll2009l ) rules out the pos- 
sibility that the majority of galaxies observed ai z — 4 have 
been forming stars at constant or declining rates since z Js 5. 
This is because constant or declining SFHs would predict 
higher SSFRs at earlier times, in conflict with their observa- 
tions. They noted that smoothly-rising SFHs could explain 
their observations since in th is case galaxies would grow in 
Mr and SFR simultaneously. IPapovich et al.l (1201 0) studied 
the evolution of the UV luminosity and stellar mass at con- 
stant number density and concluded that observations re- 
quire smoothly-rising SFHs and are consist ent with a scale- 
i nvari ant shape from z = 8 — >■ 3. Finally, iMaraston et al.l 
(2010) argued in favor of exponentially-rising over decay- 
ing SFHs as models for star-forming galaxies at z ~ 2 for 
two reasons. First, they found that rising SFHs provide bet- 
ter fits both to observations and to predictions from semi- 
analytic models of galaxy formation. Second, rising SFHs 
naturally explain higher-redshift samples as the fainter pro- 
genitors of lower-redshift samples. 

These suggestions exemplify how statistical samples at 
multiple epochs can constrain the mean SFH even if the 
SFHs of individual galaxies are unconstraine d. In this work, 
we us e this idea to build upon the results of IFinlator et al.l 
((20071) by asking whether smoothly-rising SFHs can explain 
the observed statistical properties of galaxies at z = 6-8. 
Broadly, since all of our simulated high-redshift galaxies ex- 
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msPH mdark 






12 0.469 0.235 1.19 15.0 

24 0.938 1.87 9.55 120 

48 1.88 15.0 76.4 962 



''Box length of cubic volume, in comoving 

h^^Mpc. 

''Equivalent Plummer gravitational soft- 
ening length, in comoving h~^kpc. 

■^AU masses quoted in units of IO^Mq. 

'^Minimum resolved galaxy stellar mass. 

perience smoothly-rising SFHs, any comparison with obser- 
vations can be viewed as an indirect test of these models. In 
most cases, however, the constraining power of large sam- 
ples at multiple epochs will render our comparisons direct 
tests of the smoothly-rising scenario. 

We describe our simulations in Section (2] In Section [S] 
we review the rising SFH scenario and summarize its predic- 
tions. We compare our predicted UV LFs with observations 
in Section m We perform detailed comparisons between sim- 
ulated and observed SEDs in Section [S] Additionally, we use 
our simulations to predict the completeness of current 2^6 
surveys. In Section [6] we study the predicted relationships 
between stellar mass, star formation rate, star formation 
history, and metallicity, comparing with observations where 
possible. In Section [71 we use our predicted halo occupation 
distribution to argue that recent clustering observations are 
consistent with smooth SFHs, and we discuss further tests 
of our model and alternative bursty SFH scenarios. In Sec- 
tion [S] we comment on processes that cause SFHs to depart 
from our smoothly-rising, scale-invariant scenario at lower 
redshift. Finally, we summarize our results in Sec [9] 



2 SIMULATIONS 

2.1 Numerical Simulations 

We ran our cosmological hydrodynamic simulations using 
our custom version of the parallel cosmological galaxy for- 
matio n code Gadget-2 IjSpringel fc Hernauistll2002l : ISpringell 
[2OO5I). This code uses an entropy-conservative formulation of 
smoothed particle hydrodynamics (SPH) along with a tree- 
particle-mesh algorithm for handling gravity. It accounts for 
photoionization heating starting a t z = 9 via a spatiall y 
uniform photoionizing background (IHaardt fc Madaull200l[ ) . 
Gas particles undergo radiative cooling under the assump- 
tion of ionization equilibrium, where we account for metal- 
line co oling using the coUisional io nization equilibrium ta- 
bles of ISutherland fc Dopital (|l993l ). 

The assumption of ionization equilibrium with a uni- 
form ionizing background may become increasingly inap- 
propriate as we push our predictions beyond z = &, where 
cosmol ogical reionizati on is believed to have ended (for ex- 
ample. If^ et al.ll2006l '). We argued in DFO06 that the ob- 
servable galaxies at z = 6-10 live in sufficiently biased re- 
gions that their environments were probably reionized early 



rp v, 1 o- 1 i- relative to the cosmological mean, with the implication that 

Table 1. Simulation paramc^ cro. 1 1 ■ ' ^ 

the Haar dt fc Madaul l|2001l ) ionizing background is not a 

strong approximation. In future work, we will relax this 
assumption through self-consistent radiative hydrodynamic 
simulations. However, we retain it in our present work for 
simplicity, in essence adding this to the list of approxima- 
tions that our comparisons test. 

We endow dense gas particles with a subgrid two-phase 
interstellar medium consisting of hot gas that condenses via 
a thermal instability into cold star-forming clouds, which 
are in turn evapora ted back into the hot phase by super- 
novae (JMcKee fc Os trikcr 1977]). The model requires only 



one physical parameter, the star formation times cale, which 
we tun e in such a way that it reproduces the iKennicuttI 
(|l998ah relation (jSpringel fc Hernauistll20ol ). 

We have improved our treatment of the formation and 
transport of metals. We sum marise these changes here; fo r 
details we refer the reader to lOppenheimer fc Davg l|2008h . 
Star-forming gas particles self-enrich owing to Type II su- 
pernovae as before, but we now incorpor ate the metallicity- 
dependent Type II supernova yields from lChieffi fc Limongil 
(2004 ) assuming a Chabrieg (120031 ) IMF. We account for 
mass l oss from AGB stars using the delayed feedback ta- 
bles of iBruzual fc Chariot! (|2003l ) . We account for energy 
and metal feedback from Type la supernovae (both prompt 
and delayed populations). Finally, we track the enrichment 
rates of C, O, Si, and Fe separately rather than tracking 
only the total metal mass fraction. Gas particles stochasti- 
cally spawn star particles via a Monte Carlo algorithm in 
such a way that the growth of stellar mass in star particles 
reflects the underlying star formation rate in the gas parti- 
cles. Each star particle inherits the metallicity of its parent 
gas particle. 

Galaxy formation models must invoke galactic- 
scale outflows_Jn_order to avoid the overcooling prob- 
lem (jWhite fc Frenklll99H ). We have previously shown that 
models coupling the outflow speed and mass-loading factor 
(that is, the mass of material expelled per unit mass of stars 
formed) to the halo mass are uniquely successf ul in repro- 
ducin g a wide variety of observ a tions of galaxies |Dave et al.l 



I2OO6I : iFinlator fc Pavel I2OO8I: lOppenheimer et al 



and the IGM JOppenheimer fc Dav^ |2006| . 1200^ l2009al : 



lOppenheimer et al.l [2009bj). In this work, we again use 
momentum-driven outflows with a normalization ctq = 
150 kms~^. In order to improve the fidelity with which 
we implement momentum-driven outflow scalings, we now 
compute mass loading factors using the velocity disper- 
sion of each gas particle's host galaxy rather than the lo- 
cal gravitational potential. The velocity dispersion is in 
tur n computed fro m the host galaxy's baryonic mass us- 
ing lMo et al.l (| 19981 ). For furth er details on the physics treat- 



ment s in the simul ations, see Oppenheimer fc Davg (|2006l . 



200g); Oppenheime r et al.l (|2010l ) 



Our fiducial simulation volume is a cube 24/i~'Mpc long 
on each side that uses 512"^ dark matter and star particles. 
With a mass resolution limit of 128 star particles, this im- 
plies that our fiducial volume resolves stellar populations 
more massive than 1.20 x 10* A/0. This mass limit translates 
into a limiting UV magnitude via the simulated star forma- 
tion histories (SFHs) and the assumed stellar population 
synthesis model; we will show in Section [4] that this reso- 
lution limit is well-matched to current observational limits. 
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Throughout this work, we will evaluate resolution conver- 
gence by comparing with results from additional simulations 
whose volumes span 12 and 48/i~ Mpc with the same num- 
ber of particles, for an implied stellar mass resolution limit 
of 1.50 X IO^Mq at our highest resolution (see Table [T]). 

We assume a cosmology where Q.m = 0.28, JIa = 0.72, 
Ho = 70 kms"^ Mpc"\ as = 0.82, and Q,b = 0.046. 
Note that we do not correct observed photometry to our 
assumed cosmology because the implied corrections are typ- 
ically small (~ 0.05 mag) compared to the uncertainty in- 
herent in photometric redshifts (« 0.2 mag). We measure 
simulated UV luminosities using a narrow boxcar filter at 
1350 A that is not contaminated by Lyman-a emission or 
absorption by the intergalactic medium along the line of 
sight, and we do not fc-correct observed UV luminosities to 
1350 A because these corrections are also small. 



2.2 Identifying Simulated Galaxies 

We identify galaxies within our simulations as 
gravitationally-bo und lumps of s tar and gas particles 
using SKID fsee iKeres et al.l |2005| for a description). We 
sum the SFRs over each galaxy's gas particles to obtain its 
instantaneous SFR. We obtain its star formation history 
by examining the ages of its star particles. We compute its 
stellar metallicity by performing a mass-weighted average 
over its star particles. We define its gas metallicity as the 
SFR-weighted mean metallicity over its gas particles in or- 
der to mimic the metallicities that would be measured from 
its emission lines if this were possible. In detail, our use of 
the mass- weighted stellar metallicity may not be completely 
appropriate for comparison with high-redshift observations 
given that the observed luminosities are dominated by 
the youngest stars. However, we will justify this choice in 
Figure[TT]by showing that the SFR-weighted gas metallicity, 
which dictates the metallicity of the most luminous stars 
in each galaxy, typically exceeds the mass-weighted stellar 
metallicity by only 0.2 dex. 

Note that we use SKID to identify the galaxies at each 
redshift snapshot independently. This means that our anal- 
ysis does not depend on our ability to identify each galaxy's 
progenitors and descendants (except for the illustrative evo- 
lutionary trends in Figures flOl and fTT|) . Consequently, we 
refer the "SFH" at a given epoch to the stellar age distribu- 
tion, independent of whether those stars formed in-situ or 
in smaller progenitors that were accreted at earlier times. 
While this approach pr events us from studying our sim- 
ulated merger histories, IGuo fc White! (l20_08) have shown 
that mergers are expected to be subdominan t in modulating 
the S FHs of high-redshift galaxies (see also ICattaneo et al.l 
l20ld ). The generally smooth SFHs predicted by our model 
support this view. 

We obtain each simulated galaxy's stellar con- 
tinu um by convolvin g its stellar population with 
the ISruzual fc Charloj (|2003|) stel la r po pulation syn- 
thesis models assuming a IChabrien l|2003l ) initial mass 
function (I MF) and the 'Padova 199 4' stellar evolution 
models (see iBruzual fc CharlotI 120031 for references), in- 
terpolating within the tables to the correct metallicity 
and age for each star particle. We account for du s t usin g 
the foreground screen model of ICalzetti et al.l l|2000f ). 
where each galaxy's assumed reddening E[B — V) is 



tied to t he mean stell a r me tallicity using a local cal- 
ibration (JFinlator et al.l l2006h . We treat the effect of 
absorption by the intergalactic m edium ( IGM) along the 
line of sight to each galaxy using iMadaul (1995). Finally, 
we convolve the resulting synthetic spectra with the 
appropriate photometric filter response curves to produce 
synthetic photometry. We have explicity checked that our 
photometry is unaffected when we use the updated popu- 
latio n synthesis models o f Chariot fc Bruzuaj ((20071) rather 
than iBruzual fc CharlotI (12003). This is expected because 
the observed Spitzer/IRAC [3.6] flux is not sensitive to the 
contribution of thermally pulsating asymp totic giant branch 
stars at redshifts greater than 6 (see also lStark et al.ll2009l : 
iLabbe et al.ll2010l ). Note that, throughout this work, we use 
[3.6] and [4.5] to denote Spitzer/IRAC Channels 1 and 2, 
respectively, and we use Vios, J125, and Hu,o to denote the 
WFC3/IR F105W, F125W, and F160 bands, respectively. 

In an improvement over our previous work, we now 
add nebular line emission to our simulated spectra be- 
fore a£comiting_^oi__dust and the IGM opacity follow- 
ing Schae rer fc de Barrod ([20091). We first compute the ion- 
izing continuum luminosity directly from the stellar contin- 
uum, which automatically accounts correctly for the self- 
consistently predicted stellar population age and metallic- 
ity. We then compute the luminosity in all hydrogen re- 
combinat ion transitions f rom a n upper state with n ^ 50 
using Storev fc Hummen (| 19951 ). This step involves several 
assumptions. First, we use Case B recombination theory 
and omit Lyman-a. Second, we assume an electron density 
rie = 100 cm~^ and temperature Te = lO^K. This is rep- 
resentative because the luminosities vary by ^ 30% in the 
range Tie -= 10^-10"^ cm~^ and Te =1 000-30,000K. Finally, 
we follow ISchaerer fc de Barrod (l2010|) in assuming that all 
ionizing photons are absorbed by hydrogen rather than dust. 
We th en add the dominant allowed and forbid den metal lines 
using lAnders fc Fritze-v. Alvenslebenl (2003), interpolating 
linearly to each simulated galax;y's metallicity. Incorporat- 
ing emission lines in this way allows us to make physically- 
motivated predictions regarding the impact of emission lines 
on observable galaxies at z ^ 6. 

The most uncertain aspect of our treatment for nebular 
emission involves the choice of ionizing escape fraction /esc. 
Optimally, one would like to constrai n /esc observatio nally, 
but current data do not permit this (|Ono et al.ll2010l ). We 
adopt the assumption that /esc = because it yields an 
upper bound to the effect that nebular emission can have on 
galajcy spectra. This is equivalent to assuming that current 
z ^ 6 samples did not bring about cosmological reionization. 
However, we will also argue that current samples probe only 
the brightest 1/3 of all star formation at z ^ 6 anyway, 
hence this assumption is not inconsistent with the view that 
star formation dominated the reionization photon budget. 



3 RISING STAR FORMATION HISTORIES 

The goal of the current work is to test two fundamental 
predictions of our galaxy evolution model: 



(i) SFHs at early times 
rising (|Finlator et al.ll2007l '): and 
(ii) SFH shapes are scale-invariant. 



are smoothly- 
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Figure 1. Simulated mean SFHs down to z = 5.5 in four bins 
of stellar mass. Top: Mean star formation rates in galaxies whose 
stellar masses log(Mt/M0) at 2 = 5.5 are, from bottom to top, 
8.2, 8.7, 9.2, and 9.7. Bottom: Tlie same SFHs but normalized 
to the respective mean SFRs. Our simulations predict that SFHs 
at early times are smoothly rising and vary only by a scale fac- 
tor. Here we include representative normalized mean SFHs from 
our comparison volumes from galaxies with the indicated stellar 
masses. The good agreement with the SFHs from our r24 volume 
suggests that our mean SFH is numerically converged. 



We illustrate these predictions in Figure \T\ The top panel 
shows the mean simulated SFHs in four bins of stellar mass. 
We have generated these SFHs directly from the stellar age 
distributions of our simulated galaxies at 2; = 6. These 
curves immediately demonstrate that SFHs are on average 
smoothly-rising prior to z = 6. The growth timescale to in 
an exponential fit of the form M* « e*' *° evolves rapidly 
from 50 Myr at z = 12 to 350 Myr by 2 = 6. 

The top panel strongly suggests that mean SFHs of sim- 
ulated galaxies of different masses vary only by a scale factor. 
We demonstrate this second prediction explicitly in the bot- 
tom panel. Here, we have normalized each SFH by its mean 
SFR, which is just the total stellar mass formed by that SFH 
at z = 6 divided by the age of the universe a.t z — 6. The 
normalized mean simulated SFHs are nearly identical, espe- 
cially after z — 12. This implies that galaxy ages (and hence 
colours) should vary only weakly with luminosity, and that 
rest-frame optical and UV luminosities (or, equivalently, the 
inferred stellar mass and star formation rate) should corre- 
late with a slope of unity. 

This self-similarity may be surprising for two rea- 
sons. First, star formation is expected to be quenched at 
low masses owing to photoionization feedback and at high 
masses owing to the onset of hot-mode accretion, hence self- 
similarity cannot extend to arbitrarily low and high masses. 
Our simulations are sensitive to these effects because they 
include an optically-thin ionizing background and model the 



formation of virial shocks in halos. However, photoheat- 
ing only suppres ses star formation in ha l os less massive 
than ~ 10* M© iJThoul fc Weinberd 1 19961 : lOkamoto et~aD 
2008 ) while virial shoc ks only suppress inflows in halos more 
mass i ve than IO^^Mr JKeres et"al]|2005l : iDekel fc BirnboimI 
l200a : iKeres et al.ll2009l ) . The host halos of observed galaxies 
at z ^ 6 fall largely between these scales (Figure [13)) . hence 
they do not probe the processes that cause departures from 
self-similarity. 

Second, it is surprising that galaxy age should be 
independent of mass because more massive dark matter 
halos are expected to have assem bled more recently at 
any redshift (jLacev fc Cold 119931) . This might lead one 
to expect more mas sive gala^cies to be younger. However, 
iNeistein et al.l (|200a) have shown that more massive dark 
matter halos are younger than less-massive halos only if each 
halo's age refers to the epoch at which its most massive pro- 
genitor assembles half of the current mass. By contrast, if 
one identifies the halo age with the epoch at which more 
than half of the current mass has been assembled into all 
progenitors more massive than a given threshold, then the 
trend reverses and more massive halos are older. The stellar 
age distribution relates more closely to the latter age defi- 
nition if star formation is inevitable in halos more massive 
than, for example, 10* M©, hence one might equally expect 
galaxies in more massive halos to be older. However, the as- 
sembled mass fraction in all progenitors at a given redshift 
depends more weakly on mass than th e mass fraction in th e 
main progenitor (Figures 3 and 4 of INeistein et al.ll200a l. 
This means that the actual dependence of SFH on mass is 
likely to be driven by feedback processes such as outfiows 
rather than refiecting underlying trends in halo assembly 
histories. For example, the momentum-driven wind scenario 
ties the outfiow mass-loading factor to the circ ular veloc- 
ity, which in turn varies as ff(z)^/^ iJMo et al.l 11998): this 
could delay early star formation in massive halos enough to 
remove the downsizing trend in halo formation histories. 

Numerical star formation histories suffer from resolu- 
tion limitations at early times when the number of gas par- 
ticles in the protogalaxies is small. Close inspection of the 
bottom panel of Figure [T] reveals that the normalized SFHs 
exhibit larger scatter at 2; > 12, which probably reflects 
resolution effects. For this reason, it is important to test 
how sensitive our predicted mean SFH is to the simulation's 
mass resolution. The blue dashed and red dotted curves in 
the bottom panel illustrate normalized SFHs drawn from 
our higher- and lower-resolution comparison volumes, re- 
spectively. They are in reasonable agreement with the mean 
SFHs predicted by our fiducial volume. This suggests that 
our predicted mean SFH is not strongly sensitive to numeri- 
cal resolution effects. Moreover, given that the SFHs appear 
converged at z < 12 and that most of the stars observed 
during z = 6-8 formed after z — 12, the early lack of con- 
vergence does not significantly affect our results. For refer- 
ence, the mean SFHs in the bottom panel are well-fit by the 
following polynomial (where t is the age of the Universe in 
Gyr): 



M, « 





0.2 - 3.2f + 10.5*2 



3.9r 



if t < 0.27 
otherwise 



(1) 
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Figure 2. Observed (coloured points) and simulated (black) 
rest-frame UV LFsat z=6 -8. The observe d LFs are taken from 
McLure et ai] 1I2OIOI') an d iBouwens et alj llJb oQb) at z = 8; 



Oesch et al 



"l^2010a^ and iMcLure et alJ ll2010i') at 



Bouwens et alJ l|2007l ') and iMcLure et alJ l(2009l ) at 



7; and 
6. The 

aggregate LFs are compiled from our three simulation volumes fol- 
lowring DFO06; the inset panel in the bottom-left panel shows the 
three full LFs. Error bars show the jackknife error at Myv = — 19 
at each epoch; this includes Poisson uncertainty a nd cosmic vari- 
ance. Predictions include dust extinction following lFinlator et al.l 
(1200a). The simulations and observations agree to within a factor 
of 3 at all redshifts, suggesting that our star formation and feed- 
back prescriptions adequately model galaxy evolution in the late 
reionization epoch. The lower heavy and upper light curves in the 
bottom-right panel show the simulated LF with and without dust 
reddening, respectively. The impact of dust is weak owing to low 
metallicities, and it grows progressively weaker with increasing 
redshift. 



4 LUMINOSITY FUNCTIONS 

In the top-left, top-right, and bottom-left panels of Figure[2l 
we compare the simulated rest-frame 1350 A LF with the 
observed UV LFs at redshifts 6-8. Before we discuss these 
comparisons, we note that our predicted LFs are uncertain 
owing to Poisson errors, cosmic variance, numerical resolu- 
tion effects, and the choice of IMF. We estimate the contri- 
bution of Poisson errors and cosmic variance using jackknife 
resamphng and show the resulting uncertainties at a repre- 
sentative absolute magnitude of M1350 — —19 in each panel; 
this is typically 20-30%. We estimate the impact of numer- 
ical effects by comparing the predicted LFs from our three 
volumes at 2; = 6 in the inset panel; lower-resolution volumes 
typically overproduce the LF normali zation by ~ 50%. Fi- 
nally, our choice of the lChab rier (2003) IMF yields luminosi- 
ties 50% brig hter for the sam e stellar mass than we would 
find from the lSalpeteJ j 19551 ) IMF. 

The normalizations of our predicted LFs agree with ob- 
servations to within a factor of three at all three redshifts, 
which is generally within the combined theoretical and ob- 



servational uncertainties. In particular, both the simulated 
and the predicted LFs may be regarded as brightening by 
0.6-0.8 magnitudes during z — 8 ^f 6. This is a nontrivial 
test of the smoothly-rising SFH scenario. To see this, note 
that, if the average star- forming galaxy at 2: = 6 had been 
forming stars at a smoothly constant or declining rate since 
it formed, then the LF would be constant or grow dimmer to 
lower redshifts. Instead, observations strongly suggest that 
objects at constant number density grow brighter with time, 
implying that they increase their SFRs even as they in- 
crease their stellar masses (|Papovich et al.l I2OI0I ). This is 
the smoothly-rising SFH scenario. The predicted evolution 
probes the shape of the mean SFH, and the good agreement 
with observations suggests that the simulated SFH shape is 
realistic. Note that this interpretation depends on the as- 
sumption that SFHs are smooth rather than, for example, 
stochastically b rightening into and then fading o ut of ob- 
served samples IjLee et al.l l2009l : IStark et al.l [2009|) ; we will 
argue on observational grounds that the SFHs are indeed 
smooth in Section [T] 

While the normalization and its evolution are in reason- 
able agreement with observations, our simulations predict 
rather steep faint-end slopes a ~ —2.0. This is similar to 
the slope of the halo mass function, hence our model pre- 
dicts that low-mass halos possess a fairly constant star for- 
mation efficiency (defined as the ratio of the SFR to the gas 
mass). Observations indicate a somewhat shallower slope of 
~ —1.7 at 2: = 6 (for example. iBouwens et al.ll2007l ). which 
could suggest that outflow strengths (or other feedback pro- 
cesses) scale more steeply with mass in reality than in our 
simulations. By contrast, observa tions at z =7-8 indica te a 
steeper faint-end slope of « —2 IjBouwens et al.ll2010bl ). in 
better agreement with our predictions. If this evolution is 
conflrmed, it could indicate that the feedback process that 
flattens the faint end of the LF at z ^ 6 is overwhelmed by 
the high gas inflow rate at earlier times, as expected in a 
scenario where the star formation and outflo w rates do not 
"catch up" to the inflow ra te until later times IjBouche et al.l 
I2OO9I : iPapovich et"ai]|2010l ) . 

These luminosity functions include the effects of dust, 
which we add in post-processing via a prescription that ties 
the (assumed) du st to the (predicted) metallicity following a 
local calibration (jFinlator et al.ll2006l ). We illustrate in Fig- 
ure [3] how the resulting colour excess E{B — V) varies with 
stellar mass and absolute UV luminosity (including dust) 
at z = 7. Broadly, we predict that E{B - V") =0-0.1 for 
the majority of objects at all luminosities although a sub- 
stantial fraction of objects brighter than Muv = — 19 may 
have E{B — V) =0.1-0.2. The tendency for brighter objects 
to be dustier reflects the (predicted) luminosity-metallicity 
and (assumed) metallicity-dust relations. 

The bottom-right panel of Figure [2] illustrates how this 
dust impacts the predicted UV LFs. Thick and thin curves 
show the predicted UV LFs with and without dust, respec- 
tively, at z = 6 (blue dashed) and z = 8 (solid red). Our dust 
model predicts that dust reddening is weak at early times 
owing to galaxies' low metallicities and the strong predicted 
mass-metallicity relation. Additionally, the dust column at 
a given luminosity is predicted to decrease only modestly 
with increasing redshift because the mass-metallicity and 
stellar mass-star formation rate relations evolve weakly. For 
example, at a fixed magnitude of M1350 = —19, E{B — V) 
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Figure 3. The predicted relationship at 2 = 7 between dust 
reddening and stellar mass (top) and absolute UV luminosity in- 
cluding the dust (bottom). The three loci in each panel reflect 
our three simulation volumes. Typical colour excesses of 0.05 are 
predicted at all luminosities, although a fraction of galaxies more 
luminous than Mjjv < —19 are predicted to have E{B — V) =0.1- 
0.2. 



increases from 0.04 to 0.06 during this epoch. We will return 
to the impact of dust reddening in Section [S] 

A key test of galaxy formation models involves the evo- 
lution of the star formation rate density. This can be inferred 
from the UV luminosity density {Cv), which is simply the 
integral of the UV LF. In Figure 3] we show the observed 
(coloured points) and predicted (dark solid curve) rest-frame 
Cv down to Muv = — 19 (top) and Muv = — 18 (bottom) 
at z =6-8; see the caption for details. The simulation is in 
reasonable agreement with observations to Muv <~ —19 at 
z — 6 while it slightly overproduces £^ at higher redshifts or 
when integrating to fainter limits. In both cases, the discrep- 
ancies were anticipated by Figure[2l where it is clear that the 
simulation preferentially overproduces faint {Muv > —19) 
galaxies. As in Figure[21 the observed discrepancies are com- 
parable to the uncertainty that is expected from Poisson 
errors, cosmic variance, the IMF, and numerical effects. 

The question of whether galaxies can contribute enough 
ionizing photons to reionize the Universe depends sensi- 
tively on how much star formation occurs in objects that 
are fainter than t he current observational Umit of M1350 <= 
-18 (for example. F^n fc Windhors^l2004l : lYaii et al.ll2006D . 
We now estimate this contribution in two ways. First, we 
directly measure the total predicted Li, as a function of red- 
shift within our r24 volume (still including the effects of 
dust). The dark, lower dashed curve in each panel indicates 
this total. This curve suggests that current observations do 
not detect more than Sj 25% of the total star formation 
at 2 ^ 6. In practice, this estimate suffers from numerical 
resolution uncertainty because the resolution limit of our 



Figure 4. Observed (coloured points) and simulated (curves; 
from the r24 volume) rest-frame UV luminosity density at z=6— 8. 
We obtain the observed luminosity densities by integrating pub- 
lished luminosity functions down to the indicated limits. T he ob- 
serv ed luminosity functions are from lMcLure et al.l 112010^ fblue) 
and Bouw ens et al.l ll2009lj) (green) at z = 8; from lOuchi et al.l 
(2009) fgr eeni and iMcLure et al.1 pOlOh ( green) at z = 7; and 
from .Bouwens et al.l 1I2OO7I ) (green) and iMcLure et al.l ||2009| ) 
(blue) at 2 = 6. The two observed densities at each redshift 
are offset for clarity. Solid and dashed curves indicate the sim- 
ulated luminosity density including bright galaxies (down to the 
indicated limits) and all simulated galaxies, respectively, while 
heavy red and light blue curves include and omit dust reddening, 
respectively. The simulations and observations are in reasonable 
agreement at all redshifts although the simulations slightly over- 
produce the luminosity density at low luminosites and high red- 
shifts, which we ascribe to resolution effects. Broadly, we predict 
that observations capture less than 20% of the luminosity density 
at 2 ^ 6. 



r24 volume corresponds to a luminosity only slightly fainter 
than M1350 <= —18. However, the rl2 volume indicates 
that the true LF continues smoothly to much lower lumi- 
nosities (Figure [2|. Hence we may obtain a second estimate 
for the total luminosity density as follows: If we assume that 
the luminosity density per unit halo mass varies slowly with 
halo mass, then we may compute the unresolved fraction 
by simply computing the ratio of collapsed matter in ha- 
los down to the typical host halo mass at M1350 = — 18 
(5 X 10^" M© at 2 = 6 in our simulations) versus the amount 
of collapsed matter in halos more m assive than the photo - 
heating feedback scale of 3 x 10* M© (|Okamoto et al.ll2008l ). 
Using a standard Press-Schechter calculation, this ratio is 
0.079. However, because our outflow model predicts that 
outflow strength should scale as M^ , the specific lumi- 
nosity in low-mass halos is lower, and the unresolved fraction 
should also be smaller. Correcting for this, we estimate that 
observations down to A/1350 = — 18 probe the brightest 37% 
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^ with a la scat- 
— >■ 6. This factor 



of the total luminosity density, in reasonable agreement with 
our direct numerical estimate. 

The lighter solid and dashed curves show how our re- 
sults change if we neglect dust extinction. Dust extinction 
suppresses the bright luminosity density by a factor of 2-4 
and the total luminosity density by ~ 50%, where the bright 
galaxies are more extincted because they are more metal- 
rich and hence more dusty by assumption. The widening gap 
between the dusty and dust-free luminosity densities from 
= 7 — >■ 6 tracks a moderate growth in the normalization 
of the mass-metallicity relation. 

We now pause to present our predicted conversion be- 
tween UV luminosity and star formation rate. For SFHs that 
are smooth on timescales longer than 1 Myr, this conver- 
sion is linear; tha t is, SFR = csfr * iuv IjMadau et al.ll 19981 : 
lKennicuttlll998br ). Calculating csfr requires knowledge of 
the galcixy's SFH, metallicity, and dust extinction. Within 
our simulations, the SFH and metallicity are predicted self- 
consistently, hence we may predict csfr in the absence of 
dust. Considering only galaxies with Muv ^ —18 and using 
a narrow boxcar at 1350 A to compute the UV luminosity, 
we find that csfr is independent of luminosity and equal 
to 2.0 X 10^* ergs s"^ Hz"^ (Mq yr~^)" 
ter of 30% t hroughout the i nterval z = 9 
assumes the IChabrien ((2003) IMF for consistency with our 
simulations. It may be compared directly with conversions 
that are used to derive star formation rates from observed 
UV luminosities at high redshift. 

Returning to our predicted star formation histories. Fig- 
ure [S] compares the simulated and observed stellar mass 
densities p, at z ^ 6, again integrating only over galax- 
ies down to Muv ^ —18. Comparing the heavy solid red 
curve with the coloured points reveals that the simulation 
is in surprisingly good agreement with observations at these 
epochs. The agreement in the normalizations may be fortu- 
itous. For example, our rl2 volume predicts a specific star 
formation rate 30% lower than the r24 volume in the range 
where they overlap, which implies that increasing the mass 
resolution of our r24 volume by a factor of eight would boost 
the predicted p, by roughly this factor. However, the agree- 
ment in the evolution of p* is nontrivial, and supports the 
smoothly-rising SFH scenario. To demonstrate this, we use 
a green dot-dashed curve to illustrate how p, would evolve 
from z = 6 — >■ 9 under the hypothesis of constant SFHs. In 
this trivial case, p,{Muv ^ —18) is constant and equal to 
the value at z — 6. Adopting the observed p, and simulated 
p* for galaxies with Muv ^ —18 (including dust) at z = 6, 
we find that the predicted evolution is significantly shal- 
lower than observed. The evolution is shallower because, in 
the realistic case of rising SFHs, galaxies become fainter and 
drop out of the sample at higher redshifts. This demonstra- 
tion is "conservative" in two senses: First, our simulations 
may overproduce p, at z = 6 (Figure [2]); repeating the ex- 
ercise with a lower p, would yield shallower evolution. Sec- 
ond, adopting exponentially-decaying SFHs would lead to 
increasing p« with increasing z (because galaxies near the 
limit Muv = — 18 at one epoch would fade below it at later 
epochs), in serious confiict with observations. Hence the ex- 
cellent agreement between the predicted and observed slopes 
is a nontrivial success of the smoothly-rising SFH scenario. 

As in the case of the UV luminosity density, we estimate 
observational incompleteness using our direct numerical pre- 
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Figure 5. Observed (coloured points) and simulated rest- 
frame stellar mass density p* in our r24 volum e at z =6- 
8. The blue h exagon is from IStark et al.l 1120091 ): the green 
hexa gon is fromlGonzalez et al.l 1120091 ): the magenta hexagon is 



from 
from 



Labbe et al 



Labbe et al. 



200£ ) ; and the red hexagon and square are 
1 20091 ) and JBouwens et al.ll2009bl) . Observa- 
tional points have been cor rected to account fo r galaxies brighter 
than Muv ^ —18 following lLabbe et al.l 1120091 ), and we have ad- 
ditionally subtracted 0.18 dex to convert to a IChab rier (200^) 
IMF. Thick red and thin blue solid curves show the simulated 
pt, to Muv <= —18 including and neglecting dust, respectively, 
while the dashed curve shows the full simulated p». The green 
dot-dashed curve extrapolates the observed p* at z = 6 back- 
wards assuming constant SFHs and the simulated star formation 
rate density at z = 6. The simulated normalization and evolu- 
tion are in good agreement with observations. Moreover, the pre- 
dicted evolution is in much better agreement than non-rising SFH 
scenarios. Current observations probe 20-40% of the total p, at 
z > 6. 



dictions. The dashed curve indicates the complete stellar 
mass density over all galaxies in the simulation. Comparing 
the heavy solid red and black dashed curves suggests that 
current observations do not detect more than 20-40% of the 
true p* at z Js 6. This is similar to our estimate of the de- 
tected fraction of the UV luminosity density, as expected 
given that M* ex M* . 



5 GALAXY COLORS 

5.1 Mean SEDs 

In this Section, we study the mean spectral energy distri- 
butions (SEDs) of our simulated galaxies in order to gain 
some intuition into their intrinsic properties and their im- 
plied selection function. We begin by showing in Figure [6] 
the mean SEDs of our simulated galaxies at z = 7 in three 
bins of Hieo (see caption for details). Examining the simu- 
lated intrinsic SEDs first (heavy blue curves), we see that 
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Figure 6. Predicted versus observed mean SEDs at z = 7 in three 
bins of AH = 1 centered on Hieo = 26, 27, 28. Simulated SEDs 
are plotted with solid, dotted, and dashed curves. Heavy blue and 
light red curves correspond to simulations without and with dust 
extinction, respectively. Triangles, squares, and pe ntagons indi- 
catc observations from UV-selected galaxies at 2: ~ 7 IILabbe et al.l 
L2OO9). Note that the set of model galaxies that falls within each 
bin depends on whether dust is present. This gives rise to an ap- 
parent tendency for the medium and faint bins to be brighter in 
the dusty models longward of Hiao owing to their redder SEDs. 



both the simulations and the observations indicate SEDs 
that vary weakly with H\qq, as expected in a scenario where 
the shape of a galaxy's SFH varies weakly with its lumi- 
nosity. All three bins show blue intrinsic UV continua, in 
good agreement with obs e rvations (Bouwens et al.l l2009al lbl. 



a gr - - 

l2010al : iLabbe et ahllJOOgl : iFinkelstein et al.ll2009l ). By con- 
trast, the simulated [3.6/i] fluxes are noticeably weaker than 
observed at all luminosities. Adding dust extinction reddens 
the simulated mean SEDs modestly without changing the 
level of agreement with observations. 

In order to consider the UV continua and apparent 
Balmer breaks in more detail, we compare in Figure [7] the 
simulated J125 — ^leo and Hieo — [3. 6] colours. This Fig - 
ure may be compared with Figure 1 of iLabbe et al.l l|2009f ). 
Once dust is included, the simulated UV continua (that is, 
the J125 — i?i60 colours) in the medium and bright bins are 
consisent with observations owing largely to their young ages 
(100-150 Myr), low metallicities {Z/Zq < 0.5; Figure [TT|) . 
and low predicted dust columns. 

By contrast, the predicted UV continua in the faintest 
bin are too red even if dust is omitted. This suggests that 
their predicted ages or metallicities Z are too high. It is not 
likely that this discrepancy can be resolved purely through 
younger ages. This is because the UV continua of Population 
II stars are too red even in the absenc e of nebular continuum 
emission for ages greater than 10 Myr IjBouwens et al.ll2010al . 
Figure 3), and observing populations with ages younger 



Figure 7. Observed (points with errors) versus simulated (points 
without errors) colours of galaxies at 2 = 7. Circles indicate the 
mean colours in the same three bins of i?i60 as in Figures[6l circle 
size scales linearly with i^ieo flux. The effects of emission lines 
(EL) and dust on the predicted mean colours are shown. The 
top axis converts J 1 25 — JJ 16O to UV continuum slope /3 follow- 
ing |Bouweris_et_alJ l|2010al V The predicted J125 — fl^ieo colours 
are quite blue as observed although the faintest bin may be too 
red, implying that its predicted metallicity is too high. The sim- 
ulated H\QQ — [3.6] is significantly bluer than observed. Nebular 
line emission and dust both increase the apparent Balmer break 
strength, but do not close the gap with observations. 



than 10 Myr is unlikely given that these objects have dy- 
namical times of ~ 100 Myr. We note that our predicted 
ages of » 150 Myr at z = 7 a gree with the recent simu- 
lation of ISalvaterra et al.l (|2010|), despite incorporating dif- 
ferent treatments for both outflows and metal enrichment. 
A more likely possibility is that Z suffers extra suppression 
in low-mass galaxies owing to a steeper scaling between the 
outflow mass loading factor rjw and the galaxy mass. For 
example, our simulations currently assume that rjw varies 
inversely with the velocity dispersion a. Adopting a depen- 
dence rjw oc (T~" with a > 1 would furth er suppress Z for 
faint galaxies because Z oc l./(l -I- rjw) (|Finlator fc Davg 
[2003). Such a model would also bring the faint end of the 
UV LF into better agreement with observations (F i gure [2l ). 
A flnal possibility, as pointed out bv lBouwens et al.l (l2010al ). 
is a small contribution from low-metallicity stars with a top- 
heavy IMF. Such a population would bluen the UV contin- 
uum as long as it did not give rise to significant nebular con- 
tinuum emission, which in turn would require a high ionizing 
escape fraction. It would be interesting to explore this pos- 
sibil ity within the self-consistent Population III treatment 
of S alvaterra et al.l (2010); unfortunately, these authors did 
not consider their predicted photometric colours. 

Our models qualitatively reproduce the observed ten- 
dency for more UV-luminous galaxies to display redder UV 
continua IjBouwens et al.ll2009al . l2010al ) owing entirely to the 
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simulated luminosity-metallicity and (assumed) metallicity- 
reddening relations. Hence the reported colour-magnitude 
trend could imply a modest trend of increasing dust redden- 
ing at higher UV l uminosity, as has previously been observe d 
at lower redshifts (iMeurer et al.lll999l : [shaplev et al.ll200ll ). 
Note, however, that the presence of s uch a trend in the z = 
7 da ta remains unclear at present IjSchaerer fc de Barrod 
|2010| ). 

Turning to the iifieo — [3.6] colours, we find that the sim- 
ulated galaxies generically exhibit significant Balmer breaks 
at all luminosities, in qualitative agreement with obser- 
vations. We previously predicted that significant Balmer 
br eaks are expected in Lyman break samples at z — 
4 IjFinlator et al.ll200q ). and Figure [7] now extends this pre- 
diction out into the reionization epoch. 

In detail, however, the observed breaks are roughly 0.5 
mag stronger than predicted. In combination with the blue 
UV continua, these colours pose a challenge fo r our simula- 
tions. This result is similar to the finding by iLabbe et al.l 
IJ2OO9I ) that linearly rising SFHs do not produce Balmer 
breaks stronger than Hieo — [3.6] = 0.3 mag whereas the 
observations require break strengths of 0.5 mag or stronger. 
One p ossible explanation is nebular line emission. For ex- 
ample, ISchaerer fc de Barrod (|2009l . |20IO|) have used stellar 
population synthesis models including treatments for neb- 
ular continuum and line emission to show that the strong 
apparent Balmer breaks observed at z ^ 6 can be mim- 
icked by nebular line emission from extre mely young (< 10 
Myr) , low-metallicity populations (see also lZackrisson et al.l 
120081 : lOno et al.ll2010l ). Our simulations provide physically- 
motivated priors for the typical ages and metallicities at 
these epochs, allowing us to predict the strength of the neb- 
ular emission lines. We find that line emission improves the 
agreement with the observed Hieo — [3.6] colour by 0.1 mag. 
We do not find that line emission can completely explain the 
observed -ffieo — [3.6] colours because our relatively evolved, 
enriched stellar populations do not yield sufficiently strong 
emission line equivalent widt hs (in contrast to the wid e 
range of models explored by ISchaerer fc de Barrod 12003 ) . 
Note that this conclusion is conservative in the sense that, in 
modeling emission lines, we have assumed an ionizing escape 
fraction of zero. For sufficiently large ionizing escape frac- 
tions that observed galaxies could contribute significantly to 
reionization (10-50%), the lines would be correspondingly 
weaker. 

Another possible explanation for t he observed c olours 
is inherently bursty SFHs. For example. iLabbe et al.l (|2009l ) 
show that only highly bursty SFHs can simultaneously re- 
produce the observed J125 — -ffieo and ffieo — [3.6] colours at 
low masses (and then only barely). Encouragingly, they also 
find that such a model would not overproduce the scatter 
in the observed SFR-A/* relation. Hence it is possible that 
boosting our mass resolution would resolve more minor per- 
turbations and instabilities, giving rise to more bursty SFHs 
and redder predicted -ffieo — [3.6] colours. We can explore 
this possibility by appealing to the tendency for -ffieo — [3.6] 
to vary weakly with luminosity (Figure [S| and examining 
the trend at lower masses in our higher-resolution volume. 
We find that increasing our mass resolution by a factor of 
eight boosts ffiso — [3.6] by 0.02 mag on average. This is 
quite small compared to the 0.5 mag discrepancy with ob- 
servations, hence if resolution limitations are to blame then 



overcoming them may require "zoom- in" simulations of over- 
dense regions given that our current simulations already 
probe the limits of what is computationally feasible. 

It is also possible that the discrepancy owes to a patchy 
dust scenario in which a fraction of the lines of sight through 
an LBG's interstellar medium (ISM) are optically thick to 
UV; this would suppress the rest-frame UV much more 
strongly than the optical. However, it is difficult to recon- 
cile this scenario with the blue observed UV continua, which 
seem to imply that these galaxies are essentially dust-free. 

Is it possible that a different treatment for galactic out- 
flows could suppress the simulated SSFR, thereby deepening 
the predicted Ba lmer breaks? As can be seen from Figure 
2 of iDavg l|2008h . the predicted SSFR is relatively insen- 
sitive to the details of the outflow treatment. This is be- 
cause changing the outflow strength broadly increases or de- 
creases M* and SFR together, without subst a ntially chang- 
ing t heir ratio (see also iDutton et al.l 12003 : iBouche et al.l 
I2OO3 ). Hence it is unlikely that the predicted -ffieo — [3.6] 
colours could be made redder by modifying our outflow 
treatment. 

Finally, it cannot b e ruled out t hat t he photometric 
uncertainties reported bv lLabbe et al.l (120091 ) underestimate 
the true errors because the stacked fluxes could be biased by 
the presence of bright outliers. For example, if (despite their 
careful analysis) the rest-frame optical flux of a fraction of 
the sources in each bin were contaminated by incompletely- 
subtracted neighboring objects, the resulting error in the 
stacked SED would be difficult to detect given the small 
number of sources in each bin. This reinforces the need for 
larger samples both in order to reduce errors and in order 
to allow for more accurate uncertainty estimates. 

In summary, our simulations produce blue UV continua 
and significant Balmer breaks at all luminosities, in qual- 
itative agreement with observations. The faintest observed 
galaxies are 0.1 mag bluer than predicted, which may sug- 
gest that our feedback model insufficiently suppresses the 
metallicities of low-mass galaxies. Our predicted UV con- 
tinua qualitatively reproduce the reported colour-magnitude 
trend as long as we include our metallicity-dependent dust 
reddening prescription. The simulations produce apparent 
Balmer breaks that are weaker than observed, even if we 
correct for resolution effects and account for nebular line 
emission. This is the most dramatic discrepancy that we 
have found, and it emphasizes the need for a better obser- 
vational understanding of the nature of the observed ap- 
parent Balmer breaks at z ^ 6 through A'-band imaging, 
mid-infrared spectroscopy, or mid-infrared imaging with im- 
proved spatial resolution. 



5.2 Direct Spectral Energy Distribution Fitting 



In Section 15.11 we compared the predicted and observed 
colours of galaxies at 2: ~ 7 without allowing the redshift 
or the amount of dust extinction to vary. This approach 
is justified if the photometric redshifts are accurate and 
the dust extinction is negligible. On the other hand, given 
that we do not model the dust extinction self-consistently 
and that the redshift is not well-constrained, a complemen- 
tary approach is to allow these parameters to vary freely 
and ask how well the simulations can account for observa- 
tions in principle. In this Section, we use our Bayesian SED- 
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Figure 8. Black points with errors denote the stacked SEDs of 
2 ~ 7 LBGs from lLabbe et al.l 12009 ): the -ffigo magnitude ranges 
for each bin are indicated. The blue curve indicates the spectrum 
of the best-fit numerically-simulated galaxy. Magenta crosses in- 
dicate synthetic photometry of that spectrum, offset slightly in 
wavelength for clarity. The simulated spectra yield reasonable 
agreement at all wavelengths, although the WFC3/IR measure- 
ments in the brightest bin are not well-fit. Optical emission lines 
readily explain the observed [3.6] — [4.5] colours. 



fitter SPOC (JFinlator et aLllJOOTI ) to determine how well the 
smoothly-rising SFHs fit observed SEDs at z ~ 7 and discuss 
the implied physical properties. 

Here we brie f ly re view SPOC; for details and tests 
see iFinlator et al.l (|2007l ). SPOC operates similarly to other 
SED-fitters. Starting with a set of model SFHs, it generates 
a library of SEDs by resampling the SFHs over a grid of red- 
shift z and dust column Ay- Note that, in an improvement 
over our previous work, we include optical emission lines as 
described in Section [2. 21 After comparing the library SEDs 
with the measurements, it uses a Bayesian analysis to derive 
the posterior probability for the various physical parameters 
from the likelihoods (exp(— x'^/2)) and any priors. The nov- 
elty of this approach is that we directly adopt numerically 
simulated SFHs and metallicities rather than assuming toy- 
model SFHs (for example, constant or declining). This im- 
poses physically-motivated priors on the results because the 
prior probability of a given combination of age, metallicity, 
M* , and SFR is proportional to the number of model galax- 
ies that form with these parameters. It is for this reason that, 
in contrast to conventional SED-fitting studies, SPOC only 
resamples the model SFHs in z and Ay- The number of free 
parameters is difficult to define given that the simulations 
do not sample parameter space uniformly. However, given 
that metallicity and SFR correlate tightly with M» whereas 
age varies only with redshift, there are effectively only three, 
namely M, , z, and Ay- The resulting likelihoods are equally 
as interesting as the constraints because they quantify how 



well numerical SFHs can account for the full SEDs of in- 
dividual objects, automatically identifying observed colours 
that challenge the models. 

As a test of how well our numerical simulations can 
account for the full SEDs of reionization-epoch galaxies, 
we apply SPOC t o the three binned SEDs at z ~ 7 
from lLabbe et al.l (|2009f ) in order to take advantage of their 
high signal-to-noise. We derive our model SFHs from snap- 
shots of our 24 /i~^Mpc volume at a number of discrete 
redshifts betw een z = 6-8. We do not consider low-redshift 
solutions (see iMcLure et al.ll2009l V Perturbing each model 
galaxy in redshift and sampling 100 different values of Ay 
between 0-1.5, we generate a library of 5 million model 
SEDs for comparison. We treat detections with less than 
2(T significance as 2a upper limits and adopt the reported 
measurements and errors in the other bands. We adopt the 
WFC3/IR Yio5 and VLT/ISAAC Ks bands to compute sim- 
ulated Y and K fluxes, respectively. 

In Figure O we compare the observations (black points) 
with the spectrum of the best-fit combination of model SFH, 
redshift, and Ay (blue curve) and its associated photome- 
try (magenta crosses) . Examining the brightest bin first (top 
panel), we find that our best-fit SED matches the observed 
fluxes in all bands other than J125 and -ffieo, both of which 
it underproduces by ~ 2a. Intriguingly, we find that opti- 
cal emission lines completely explain the observed K — [3.6] 
and [3.6] — [4.5] colours. Many reionization-epoch SEDs show 
evidence for a bl ue optical continuum suggestive of opt i- 
cal emission lines r Labbe et al.ll2009l : iGonzalez et al.ll2009l ). 
Our modeling suggests that these colours are real rather 
than reflecting observational uncertainty (for example, in 
the IRAC measurements), with several implications. First, 
it reinforces the ne ed for SED-fitting studies to accou nt for 
nebular emission (jSchaerer fc de Barrod |2009| . 120101 ) . Sec- 
ond, it surprisingly suggests that the difficulty with the ob- 
served Hiao — [3.6] colours (Figure [7}, if observational in 
origin, may owe as much to the WFC3/IR as to the IRAC 
bands given that the K — [3.6] and [3.6] — [4.5] colours are 
not difficult to reproduce. Finally, it reinforces the need for 
deep K-hsmd imaging in order to constrain the strength of 
the Balmer Break. 

We may derive constraints from our posterior probabil- 
ity distributions by marginalizing over all but one parame- 
ter at a time and determining the resulting 68% confidence 
intervals. This method generally suffers from significant de- 
generacies between the derived age, metallicity, and dust 
extinction (for example. IShaplev et al.ll200ll : |Papovich et al.l 
,20Ql!). However, restricting our attention to the SFHs that 
arise within our simulations partially breaks these degen- 
eracies because the predicted range of SFHs i s narrow, effec- 
tively introducing tight priors (for details, see IFinlator et al.l 
l2007f ). We find 2 = 6.7±0.1, Ay = 0.6±0.1, log(M./M0) = 
9.6 ±0.1, M. = 16.O±3.8M0yr-\ and Z = 0. 007 ±0.0009 
The in ferred SFR is 60% larger than inferred bv lLabbe et al.l 
120091), and the inferred stellar mass is 10% lower (both num- 
bers are corrected for the different assumed IMFs). Mean- 
while, the inferred dust column is much larger than would 
be expected from the blue J125 — Hmo colour. This explains 
why the model SED is too red in Hiao — [3.6] whereas, with 
our fiducial dust prescription, it was too blue (Figure[7|. Not 
surprisingly, the associated J125 — -ffieo colour is now also 
too red whereas it was in good agreement with the observed 
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colour in Figure [T] These constraints are driven largely by 
the A'-[4.5] bands, and the lower stellar mass is permitted in 
part by the in clusion o f nebular emission lin es (for example, 
ISchaerer fc de Barros 201ol : lOnoet alJlJOlol ). 

In the bottom panels, we compar e the stacked SEDs 
from the medium and faint bins of iLabbe et al.1 (|2009l ) 
with the best-fit model SEDs. These bins yield generally 
better fits to our models than the brighter bin. As be- 
fore, the simulated SEDs readily reproduce the observed 
[3.6] — [4.5] colours owing to the presence of optical emis- 
sion lines. The medium stack implies the physical parame- 
ters z = 6.8±0.1, Av = 0.08±0.08, log(M./MQ) = 8.8±0.1, 
SFR=3.1±O.4M0yr-\ and Z = 0.0040 ± 0.0006, while the 
faint stack corresponds to z — 6.9 ± 0.1, Av ~ 0.09 ± 0.08, 
log(M./M0) = 8.5 ± 0.1, A/. = 1.5 ± O.3M0yr~\ and 
Z = 0.0032 ± 0.0007. As before, the stellar masses are 
lower than would be inferred without emission lines. Mean- 
wh ile, the infer r ed SF Rs are 10-30% lower than inferred 
by iLabbe et al.l (|2009|) (after accounting for the different 
IMFs) owing to the low predicted metallicities. The inferred 
dust columns are significantly lower than in the bright bin. 
Note, however, that deeper observations in K may eventu- 
ally demand higher dust columns in the faint bins as well. 



5.3 Survey Completeness 

We may use our simulated SEDs to predict the completeness 
of current surveys. In detail, of course, completeness varies 
with the selection strategy. However, a tight correlation be- 
tween stellar mass and SFR implies that the completeness 
of UV-selected samples is dominated by the det ection limit 
rather than by colour cuts (jFinlator et al.ll2006r ). This is es- 
pecially true when the dust column is low. For example, the 
colour cuts in lOesch et al.l (|2010a|) do not eliminate any of 
our simulated galaxies at z = 7. 

In Figure|9] we plot using heavy curves the number frac- 
tion of galaxies that arc selected by the 5cr detection limits in 
Ji25 and HiQo of Bouwons ct al. (2010ai) versus stellar mass 
both with and without dust. Galaxies with M* ^ 10^ M© 
are selected efficiently assuming that dust reddening is small 
{E{B ~V) ^ 0.05) owing to their high star formation rates. 
Dust increases the 50% completeness limit from IO^'^A/q 
to lO^'^Af© by extincting galaxies whose UV luminosities 
are already near the detection limits. Galaxies with stellar 
masses less than IO^'^A^q are not selected efficiently regard- 
less of dust because their star formation rates are too low. 
The 5a detection limit for a 10^ second exposure through the 
JWST/NIRCAM F200W fiher will be AT = 29.72; Q we show 
how this improves the mass completeness using light curves. 
Clearly, a relatively shallow JWST exposure will readily de- 
tect objects that are 10% as massive as the least massive 
objects detected in the WFC3/IR ERS fields. The impact 
of dust extinction will be weak owing to the low metallicity 
(Figure [TT|) . 

These results are robust to our choice of input physics 
or mass resolution. This is because the mass completeness 
is mostly sensitive to the specific star formation rate, which 
in turn is relatively insensitive both to our choice of outflow 
strength and mass resolution. 
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Figure 9. Completeness versus stellar mass at z = 7. Curves 
show how the number fraction of galaxies that are selected by 
two different criteria varies with stellar mass in simulated catalogs 
with (solid red) and without (dashed blue) dust reddening. Heavy 
curves correspond to the current 5cr limiting m agnitude in the 
WFC3 /IR ERS fields of J ^ 27.7 A H s^ 27 A JBouwens et al.l 
l2010al ). Light curves correspond to the 5cr detection limit for a 
10,000 second exposure in the JWST/NIRCAM F200W band. 
Current surveys detect most galaxies more massive than IO^A/q 
and few that are below 10^-^ Mq, while JWST will easily detect 
objects that are 10% as massive. 



6 PHYSICAL PROPERTIES 

In this section, we summarize the predicted physical proper- 
ties of our simulated galaxies at 2 ^ 6, where possible com- 
paring with observational inferences. The most important 
physical process impacting these predictions is our simula- 
tion's implementation of momentum-drive outflows, and we 
have previously explored the effects at 2: = 6 of varying our 
outflow treatment in DFO06. 



6.1 Star Formation Rates 

We show in Figure [10] how SFR relates to A/* in our sim- 
ulated catalogs (black points) at four different redshifts. 
Broadly, we predict SFRoc A/P'^^^^'^^ at aU redshifts, con- 
flrming our previous results (DFO06). The predicted slope 
at 2 = 7 of 1.05 is in excellent agreement with the reported 
slope of 1.06 ±0.1 (Labbc ct al. 2009), and the agreement is 
good at the other epochs as well. This comparison is a non- 
trivial test of our model, and given that our simulations were 
not tuned to reproduce s uch a trend — in fact, this class of 
simulations predicted it (jDave et al.l I2OO0I : [Weinberg et al.l 
2002) — the level of agreement is remarkable. 

There are two possible interpretations for a near-unity 
slope. One possibility is that galaxies begin forming stars at 
different epochs, but once they begin, their growth is expo- 
indaae-ftiiai with a timescale A^*/Af*. This timescale is observed 
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Figure 10. Panels a— c show simulated (black points) and ob- 
served (coloured shapes) SFR-M, relationships. The three simu- 
lated loci correspond to our three volumes, and their good agree- 
ment suggests that numerical errors are weak. We include the 
slope and dispersion about the best-fit linear relation from the 
r24 volume at each redshift. Panel d shows the simulated rela- 
tionships at ^ = 9 (red crosses) and z = 6 (blue points) as well 
as typical evolution of individual galaxies from 2 = 9 — » 6 (ar - 
rows). The red square at 2 = 8 is from iBradlev et al.l ll200a ). 
where we assume constant star formation and Z« = 0.0004; the 
black tri angle corresponds to the best-fit SFH to the 2 = 7.7 
stack of iLabbe e t al.l 112 0091): a nd the magenta hexagon is Ob- 
ject 1 from ICapak et al.l ||2009|V At 2 = 7, black triangles are 
the s tacked result s from iLabbe et al.l 1 120091) ; the red square is 
from lEgami et al.l l|200a ) , where the point represents the aver- 
age of the inferences assuming a 100 Myr exponentially decay- 
ing SFH over the different assumed metallicities, and the er- 
rors reflect this error added in quadrature to assumed photo- 
metric errors of (5%, 20% ) in the (SFR, mas s); the magenta 
hexagon is Object 3 from ICapak et al.l 1120091); and the green 
square is the 2 = 6.6 LAE stack of lOno et al.l l|2010r ). where we 
plot result from the "pure stellar" model with Z/Z p, = .2. At 
2 = 6, the blue pentagon is LA34 from iLai et al.l (|2007J); the 
red square is from iDo w-Hvg elund et al.l 1120051 ); the black trian- 
gles are from I E vies et al.l (|2007|); and the gree n square is from 
the stacked SED at 2 = 5.7 of lOno et al.l 1 120101 ) assuming "pure 
stellar" models with Z/Zq = 0.2. We have subtracted 0.18 dex 
from all observationally inferred SFR and M* to adjust from the 
Salpeter to the Chabrier IMF. Our simulations predict that SFR 
oc Mt, with tight scatter and a weakly-evolving normalization. 
The slope and evolution of the simulated trend are in good agree- 
ment with observations, but the normalization may be offset by 
a factor of 2—3. 



to be ^ 500 Myr at z = 6 (jStark et al.ll2009l ), although our 
simulations indicate somewhat shorter timescales of 200- 
350 Myr (note that accounting for optical emission lines 
will bring observations into closer agreement with our pre- 
dictions; Section 15. 2p . The short growth timescales could 
then imply a bursty star formation scenario in which galax- 
ies "turn on" at different times, form stars with a short 



duty cycle and a constant exponen tial growth timesc ale, and 
then fade into a quiescent state (jStark et al.ll2009l ). There 
are two difficulties with this picture. One is that it is not 
obvious why galaxies that begin forming at very different 
times should nonetheless obey the same exponential growth 
timescale with negligible scatter. The other is that bursty 
models predict a large populatio n of quiescent galaxies at 
z_^4 that has not been observed (|Brammer fc van DokkumI 
[20071). We will return to these points in Section [T] 

A second interpretation of the near-unity slope is that 
SFHs at z > 6 have a scale-invariant shape. This is be- 
cause, if all galaxies begin forming stars at the same time 
and with the same SFH shape, then the SFR and the M» 
both vary linearly with an overall scale factor. Consequently, 
the slope remains near unity until a mass-dependent pro- 
cess such as hot-mode accretion or AGN feedback breaks 
the scale-invariance and flattens the slope (Dave 2008). In 
this smoothly-rising SFH scenario, the growth timescale is 
dominated by smooth inflows, which in turn are regulated 
by the competition between the growth of halo potential 
wells and the decre ase in the cosmic density iJBouche et al.l 
1 20091 : iDutton et al.l i2009). It is not easy to distinguish ob- 
servationally between the bursty and the smoothly-rising 
scenarios using UV-selected samples, although our models 
(and, for that matter, all hydrodynamical simulations) sup- 
port the latter view. However, we will discuss other obser- 
vational probes in Section [T] 

The predicted normalization may be slightly offset from 
observations. This offset could have two possible implica- 
tions. First, it suggests that observationally inferred stellar 
masses are too high because they do not include the effect 
of optical emission lines. Accounting for this would likely re- 
duce observationally inferred stellar masses by 0.1-0.3 dex 
(Figures [7] and [8l , improving the agreement with predic- 
tions. In support of this view, Labbo ct al. (2009) have found 
that applying an observationally-motivated estimate for the 
strength of the OH emission line lowers the typical inferred 
Af, at z = 7 by 0.17 dex. The second implication is that our 
simulated SFRs are slightly too high, which we previously 
noticed evidence for in Figure [2] It is tempting to suppose 
that strengthening our outflows would solve this problem. 
Unfortunately, as noted in Section [5l this would not sup- 
press the predicted SSFR because increasing the outflow 
mass-loading factors reduces M« and M« togeth er without 
changing their ratio (see also lDutton et al.ll2009l ). Hence it 
is not simple to understand how the bulk of the offset could 
be attributed to our input physics. 

The predicted dispersion about the mean trend is 
« 0.1 dex at all redshifts. This tight scatter reflects the 
tendency for star formation to be driv en by smooth gas 
accretion rather tha n by interactions (JKeres et al.l l2005l : 
iBirnboim et al.ll2007l ). It is some what less th an the reported 
scatter of 0.25-0.3 dex at z ~ 7 (JLabbe et a l. 2009). Given 
that observational uncertainties boost the observed scatter, 
it is possible that the true scatter is consistent with our 
predictions. It is also possible that our limited numerical 
resolution does not account fully for the minor interactions 
or discrete infalling clouds that perturb galaxies away from 
their equilibrium SFR. We evaluate this possibility by plot- 
ting the predicted SFR-M* trend from all three of our sim- 
ulation volumes at each redshift (hence the three separate 
"clumps" in the simulated locus) . These three volumes span 
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a factor of 64 in mass resolution and should immediately 
reveal which predictions are strongly sensitive to numeri- 
cal effects. The offsets in the trends are small compared to 
the scatter, and the scatter does not vary with resolution. 
This indicates that the predicted trends are numerically con- 
verged. Moreover, it is interesting to note that the mass res- 
olution of our 12/i~^Mpc simulation is 50% higher than the 
gas-p hase mass resolution with which iMihos &: Hernguistl 
IJ199Q ) demonstrated that mergers can boost SFRs of MW- 
scale galaxies by factors of 10-100. The fact that our pre- 
dicted SFR-M, scatter remains small despite our resolving 
power supports the view that the predicted scatter is robust 
to resolution limitations, and that dramatic merger-driven 
star bursts are ind e ed un common at high redshift. 

iDutton et al.l (|2009l ') use a semi-analytic model for the 
growth of disk galaxies within a smoothly-accreting halo to 
predict a scatter of ~ 0.1 dex at all redshifts. They attribute 
this to scatter in the mass accretion histories of the host ha- 
los. They then speculate that the larger observed scatter 
could owe to scatter in the gas accretion rate at fixed halo 
mass, which their model does not account for. Our simula- 
tions account for this effect self-consistently, hence the fact 
that our predicted scatter agrees with theirs argues that the 
observed scatter is not dominated by scatter in the gas ac- 
cretion rate at fixed halo mass. 

In the bottom-right panel, we show how galaxies' 
masses and SFRs evolve from z = 9 — >■ 6 in each of our 
simulation volumes. Red crosses and blue points show the 
simulated loci at z = 9 and 6, respectively. Comparing these 
loci reveals that the SFR at a given M« (that is, the nor- 
malization of the M,-Mt relation) is predicted to decline 
by ~ 0.3 dex from 2: = 9 — >■ 6. As we have already seen, 
observations also indicate that this normalization evolves 
weakly out to at least z = 8, qualitatively supporting this 
picture l|Gonzalez et all [200 9: Stark ct al. 2009). This slow 
evolution indicates that galaxies grow roughly exponentially 
during this epoch, as can also be seen from the shallow evo- 
lution of the slope of the mean SFH after z = 9 in Figure [l] 
As pointed out bv lStark et al.l (|2009l ). the only smooth SFH 
that is consistent with these trends is a rising one. 

In order to illustrate how this works, we have added 
arrows indicating the typical evolution of individual galax- 
ies during this epoch. We trace this evolution by looping 
through the 10 most massive galaxies identified in each sim- 
ulation at 2 = 6 and searching for their most massive pro- 
genitors at 2 = 7-9. We then average the masses and SFRs 
of these 10 galaxies at 2 = 6 and their most massive progen- 
itors at 2 = 9; the resulting average evolution is represen- 
tative of how all galaxies are predicted to evolve. Galaxies 
evolve in a direction that is slightly shallower than the mean 
trend at a given redshift. The tendency to evolve nearly par- 
allel to the observed trend indicates smoothly-rising, nearly 
exponential growth while the slightly shallower slope reflects 
a slow decl ine in gas accretion rates owing to cosmological 
expansion (JBouche et al.ll2009r ). 

Observational inferences suffer from two systematic 
biases, either of which could introduce artificial offsets 
or scatter at each redshift. First, the relevant data were 
measured by many different groups (see the caption), 
each of which introduces slightly different photometric 
techniques and selec tion biases. The reported e rrors may 
be underestimated (jSchaerer fc de Barroj 120101 ) although 



side-b y-side comparisons sug gest that this problem is not 
large (JFinkelstein et al.ll2009n . Second, different groups used 



different assumptions in SED-fitting. This introduces scatter 
because the choice of SFH and stellar metallicity biases the 
inferred stellar mass by 10-30% on average, with biases of up 
to a factor of 10 possible for certain model SFHs (for exam- 



ple. IShaplev et al.ll200ll:|Papovich et al.ll200ll : lFinlator et all 



I2OO7I : iMaraston et a l. 2010). It would be of interest to repro- 
duce these observational constraints using unified modeling 
assumptions at al l epochs in order to minim ize these ef- 
fects (for example. ISchaerer fc de Barrosll2010h . Viewed dif- 
ferently, however, it is intriguing that, despite the variety of 
observational and modeling techniques that underlie these 
observations, the resulting constraints tell a consistent story: 
The observed SFR-M* trend obeys a near-unity slope, its 
normalization evolves slowly, and it may be slightly offset to 
lower SFR or higher M« than predicted by our model. It is 
not unreasonable to suppose that applying a consistent set 
of modeling assumptions to these data would only tighten 
and reinforce the inferred trends. 

In summary, current observations indicate that our sim- 
ulations reproduce the observed slope and evolution of the 
SFR-i\f, trend at 2 ^ 6, while the normalization may be 
somewhat ofi'set and the scatter is 50% lower than ob- 
served. The near-unity slope supports the view that ob- 
served reionization-epoch galaxies began forming stars at 
similar epochs and possess SFHs that difi^er from one an- 
other only by a scale factor. The offset in the normalization, 
if confirmed, indicates that the inferred stellar masses are 
too high because they do not correct for optical emission 
lines. The weak observed evolution in the normalization is 
inconsistent with any smooth SFH other than a rising one. 
The small predicted and observed scatters are broadly con- 
sistent with the view that high-redshift star formation is 
driven predominantly by gas infiows rather than mergers, 
while the fact that the observed scatter is larger may reflect 
observational uncertainties. 

6.2 Metallicities 

Direct observations of gas-phase abundances now indicate 
that galaxies ex hibit progressiv e ly lower metallicitie s at 
higher redshifts (JErb et al.|[2006l : iMaiolino et"aLl 120081 '). In 
DFO06, we showed that this occurs naturally in the hierar- 
chical growth scenario and predicted that observable galax- 
ies at 2 ~ 6 should exhibit metallicities less than one tenth 
solar. Since then, a number of works have lent qualitative 
support to this prediction by noting that subsolar metallic- 
ities yield better fits to obse rved SEDs at 2 ^ 6 than so- 
lar metallicities (for examp le, Evles et al. 20071 : IStark et al.l 
I2OO9I : iBouwens et"aLll2010al : iLabbe et al.ll2010l ). 

On the other hand, an important result from DFO06 
was the prediction that very little of the observable star 
formation at 2 ~ 6 should result in primordial-metallicity 
"Population HI" stars because metal enrichment occurs so 
quickly. Given that our updated simulations incorporate sig- 
nificantly more realistic treatments for metal enrichment 
(Section [2} , it is of interest to determine whether these re- 
sults have changed. To this end, we use this section to update 
our predicted mass-metallicity relation at 2 ^ 6 and show 
that improving our treatments for metal enrichment revises 
our predicted metallicities up rather than down, supporting 
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Figure 11. Simulated SFR-weighted gas metallicity (blue top lo- 
cus) and mass-weighted mean stellar metallicity (red bottom lo- 
cus) versus stellar mass at 2 = 7. We assume a solar p hotospheric 
oxygen mass fraction of 0.00574 llAspIund et al.ll2009l) . The three 
"clumps" indicate our three simulation volumes. The arrows in- 
dicate the mean evolution in stellar mass and stellar metallicity 
of the 10 most massive galaxies in each simulation volume from 
z = 9 —^ G. Observable galaxies are already enriched above O.IZq 
by 2 = 7. Their evolution is slightly steeper than the mean trend 
at a given epoch owing to weak evolution in the normalization of 
the mass-metallicity relation. 



our view that reionization-epoch galaxies should not exhibit 
a significant contribution from Population III stars. 

In Figure llll we show with red points (bottom lo- 
cus) the mass-weighted mean stellar oxygen metallicities 
of our simulated galaxies at 2: = 7, normalized to the so- 
lar photospheric oxygen mass fraction. Broadly, we predict 
that observed reionization-epoch galaxies possess metallici- 
ties greater than O.IZq. We also predict a tight trend with 
« 0.1 dex of scatter in the residuals, simil ar to our previ- 
ous finding at z = 2 (|Finlator fc Davell200J ). Note that the 
predicted relation for the iron mass fraction is similar but 
shifted down by 0.2 dex, reflecting the weak contribution of 
prompt Type la supernovae at early times. 

The upper, blue locus shows the SFR-weighted mean 
gas-phase oxygen metallicities for the same galaxies. These 
points predict that gas metallicities follow the same trends as 
the stellar metallicities but boosted by ~ 0.2 dex becase stel- 
lar metallicities reflect the lower metallicities that character- 
ized each galaxy's progenitors. This figure may be compared 
to Figure 5 of DFO06, where we used a similar set of simu- 
lations to predict that the gas-phase metallicities at z —6-8 
were generally larger than 1/30^©. Our current simulations 
track metal enrichment significantly more realistically, as 
summarized in Section [21 Given this abundant increase in 
realism, the fact that our newer simulations predict higher 
metallicities than our previous work underscores the conclu- 



sion that observable stellar populations at z =6-9 should not 
contain a significant stellar mass f raction below the Popula - 
tion III threshold of [O/H] < -3 (|Bromm fc Larsonl l2004l). 

We use arrows to indicate how galaxies evolve in stellar 
mass and stellar metallicity from 2 = 9 — >■ 6 in each of our 
three simulation volumes (the gas-phase metallicity is simi- 
lar, but shifted up by 0.2 dex). We compute this mean evo- 
lution using the 10 most massive galaxies in each volume at 
z = 6 as described in Section FS.ll Galaxies evolve along a di- 
rection that is slightly steeper than the me an trend at z = 7. 
This can be understood as follows (see also lFinlator fc Davg 
[2OO3): In the presence of strong outflows, gala:xy metallici- 
ties closely track the equilibrium metallicity Z = y/(l+77w), 
where y is the metal yield and rjw is the mass-loading factor, 
or the rate at which material enters the outflow divided by 
the SFR. This equilibrium metallicity owes to strong cou- 
pling between the inflow, star formation, and outflow rates. 
In the momentum-driven wind scenario, r]w shrinks with in- 
creasing mass, hence the equilibrium metallicity grows with 
increasing mass. The tendency for galaxies to evolve slightly 
more steeply than the mean trend at a given epoch reflects 
the prediction that the normalization of the mass-metallicity 
relation grows by 0.2 dex from 2: = 9 — > 6. 

To understand the increase in our predicted metallici- 
ties with respect to DFO06, note that, at a fiducial stellar 
mass of 10^ M0 and redshift z — 6, DFO06 predicted a gas- 
phase metal mass fraction of 0.002, whereas our current sim- 
ulations predict a mass fraction of 0.006 overall, with 0.005 
of this in oxygen alone. This increase owes to the various 
improvements within both our simulations and our anal- 
ysis: The newer code predicts noticeable a enhancements 
at z = 6; the adopted mass- loading factors were twice as 
large at a given ha lo mass in our p revious simulations; and 
our adoption of a IChabrien (|2003i ) IMF when computing 
the yields significantly boosts the enrichment rates. Each of 
these effects boosts the predicted gas-phase metallicities by 
0.1-0.3 dex, hence an overall factor of three increase is not 
surprising. This difference may be regarded as an estimate 
in the uncertainty owing to our imperfect understanding of 
how metals are created and transported (although we believe 
that our current simulations are the more realistic). For ex- 
ample, if we were to re-run our simulations with higher mass 
loading factors in order to improve the agreement between 
the simulated and observed UV LFs (Figure [2]), our pre- 
dicted metallicities would shrink by 0.1-0.3 dex. Broadly, 
however, this uncertainty is far too small to change the 
fundamental conclusion from Figure 5 of DFO06: Observ- 
able stellar populations at z =6-9 are robustly predicted to 
have metallicities that are well in excess of O.OIZq, hence it 
should be possible to understand their SEDs without refer- 
ence to the predicted SEDs of Population III stars. 

Despite the fact that star formation at z ^ 6 is dom- 
inated by metal-enriched gas, a residual fraction of metal- 
poor star formation is predicted to persist even in observ- 
able galaxies. To estimate this, we have computed for each 
simulated galaxy the fraction of star formation occurring in 
gas whose oxygen mass fraction falls below 10~^Zq. In Fig- 
ure [121 we show how this fraction varies with absolute mag- 
nitude at z = 6 (points) and z = 9 (red crosses). At both 
redshifts, we predict that roughly 10~^ of all star formation 
occurs in metal-poor gas, fairly independent of luminosity. 
Hence it is entirely possible that the observed reionization- 
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Figure 12. The fraction of star formation occurring in gas with 
oxygen mass fraction less than 10"'^ Zq as a function of absolute 
magnitude at z = 9 (red crosses) and z = 6 (black points). A 
residual fraction of Population III star formation is predicted to 
persist even in bright galaxies at 2 = 6. 



epoch samples contain a tiny fraction of Population III star 
formation although the impact on the total SEDs is neghgi- 
ble. In detail, comparison of the three simulated loci suggests 
that our simulations do not completely resolve the metal- 
mixing processes that determine the Population III fraction. 
Extrapolating the trend from the highest-resolution simu- 
lation, we find that galaxies with Muv = — 17 may have 
a Population III fraction as high as 1%. This population 
may eventually be detectable through its enhanced rate of 
pair production SNe or GRBs. Note that these results are 
in qualitative agreement with previous suggestions that a 
residual le vel of Population III star formation persists af- 
ter z = 6 jTornatore et al.|[2007l : IJohnsonllJoTol : iMaio et al.1 



l2010l;lTrentiet al.ll2009l } 



ISalvaterra et al.l ( 20ld ) have recently used observed 



samples a,t z = 5-10 to test the predictions of a cosmologi- 
cal hydrodynamic simulation that differs from ours in three 
major respects: First, they incorporate an explicit treat- 
ment for the transition from Population III to Population 
II star formation whereas our simulations do not treat Pop- 
ulation III. Second, they treat the formation, destruction, 
and dispersal of dust grains in detail in order to predict 
the dust reddening whereas our model ties the normaliza- 
tion of a foreground dust screen to local observations. Fi- 
nally, their simulations model galactic outflows under the 
assumption t hat the m ass-loading factor and wind speeds 
do not vary (jSpringel fc Hcrnquist 200 3) while our simula - 
tions assume momentum-driven winds (jMurrav et al.ll2005h . 
Despite these differences, their finding that the impact of 
Population III stars on observable z = 6-8 galaxies is slight 
(^ 1% of the UV luminosity for Aduv < —18) is qualitatively 



Figure 13. top The simulated LBG HOD at z = 5.5 for sam- 
ples brighter than M1350 ^ —18 (solid blue) and -19 (dashed 
red). The arrows indicate the hypothetical minimum halo mass 
for each sample (see text). The simulated HOD lies well below 
unity for low-mass halos log(M(,/MQ) < 11 owing to strong out- 
flows, and it climbs gradually owing to significant scatter in the 
relationship between halo mass and M1350. bottom The cumula- 
tive occupancy for samples satisfying the same luminosity cuts as 
in the top panel. Thick and thin curves show the cumulative halo 
occupancy with and without additional galaxies in halos that host 
more than one galaxy satisfying the luminosity cut. The predicted 
occ upancies lie with in the observationally inferred range of 0.15- 
0.6 IILee et al.ll2009h . but they reflect a suppressed HOD owing to 
outflows rather than a star formation duty cycle. 



consistent with our own results, and justifies our decision to 
neglect the associated processes in our simulations. 



7 INTERPRETING THE OBSERVED HALO 
OCCUPANCY 

Up until this point, we have used a variety of arguments 
to demonstrate that, if Lyman Break galaxy (LBG) SFHs 
are smooth, then observations indicate that they must be 
rising rather than constant or declining. We have not yet 
considered th e possibility that LBG SFHs are highly bursty. 
For example, iLee et al.l l|2009l ) recently used accurate clus- 
tering observations ai z — 4-6 to show that the fraction of 
massive halos that host LBGs (hereafter, the "occupancy") 
is likely less than unity, and interpreted this as evidence 
that L BGs have sh o rt sta r formation duty cycles. Subse- 
quently, |Starl^et_alJ (|2009l ) built upon this idea to speculate 
that the progenitors of LBGs observed at one epoch may 
not be observable at earlier epochs since the progenitors' 
SSFRs would violate the observed non-evol ving SSFR. This 
interpretation echoes iFerguson et al.l (|2002l ) , who used rest- 
frame UV-optical measurements combined with toy-model 
SFHs to demonstrate that the progenitors of LBGs at 2 = 3 
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cannot dominate the observed star formation rate density 
at z = 4 if SFHs are smooth. In this Section, we explore 
how our simulations populate dark matter halos with LBGs 
and argue that outflows naturally give rise to the observed 
occupancy, with the implication that there is no conflict be- 
tween clustering observations and the smoothly-rising SFH 
scenario. We defer a more direct comparison of our predicted 
clustering properties with observations to future work. 



7.1 Halo Occupation Distribution 

For reference, we begin by reviewing how an observed sam- 
ple's occupancy may be inferred from its LF and cluster- 
ing behavior. The connection is made by modeling the un- 
derlyin g galaxy population's halo occupation distribu tion 
fHOD: iBerlind fc Weinberg |2002| : iBullock et ai]|2002l ). or 
the mean number of observable galaxies as a function of 
halo mass. This in turn derives from the survey selection 
probabili ty, the LF of isolated halos, and t he subhalo mass 
function feng et al.ll2003l : iLee et al.ll2009l '). To be explicit, 
let us define the isolated halo LF dNg/dL{M) as the intrin- 
sic number of galaxies with luminosity between L and L + dL 
in a halo of mass M; the selection probability for galaxies 
of luminosity L as p{L); and the number of subhalos with 
mass between m and m+dm contained within a halo of mass 
M as dNsh/dm{M). Then the mean number of observable 
galaxies as a function of halo mass {Ng)M is (dropping the 
L and M dependencies for clarity): 



r°° dN f'°° 



pM 



dNsi, dNg 



The first and second terms in Equation [2] describe the con- 
tributions of central and satellite galaxies, respectively. Gen- 
erally, one assumes that p{L) can be modeled and that the 
halo and subhalo mass functions are known. In this case, 
HOD modeling consists of devising a functional form for 
the isolated halo LF and then constraining its parameters 
to reproduce the observed luminosity function and cluster- 
ing properties. For reference, we refer to the intrinsic LF of 
individual halos including the subhalo contribution as the 
conditional luminosity function (CLF). 

After constraining the HOD, it is straightforward to in- 
quire what fraction of the halos that are eligible to host an 
observable galaxy actually do so. If it is assumed that all ha- 
los more massive than the lowest-mass halo that contributes 
to the sample Mmin are eligible, then this occupancy is given 
by 






(3) 



where dNh/dM is the (parent) halo mass function. Note 
that, throughout this discussion, we use "occupancy" to de- 
note the mean number of observable galaxies hosted by halos 
within a certain mass range, including the contr ibution of 
their subhalos. By contrast, the "duty cycle" of iLee et al.l 
l|2009f) refers to the integral of the LF over luminosity for 
a single halo, neglecting subhalos. For sufficiently small as- 
sumed scatter in the luminosity-halo mass relation, these 
quantities are the same. 

Observationally, the mean occupancy turns out to be 



robustly less than unity at z =4-5 (|Lee et al.ll2009r ). Unfor- 
tunately, the physical implication of this result is unclear 
because the nature of the galaxies that inhabit "dark" ha- 
los is unconstrained. One possibility is that all halos above 
A/min do, in fact, host a star-forming galaxy and that a frac- 
tion 1 — Ng/Nh are in a temporarily quiescent state. In this 
case, the oc cupancy Ng/Nh constrains the star formation 
duty cycle JGiayalisco fc Dickinsoiil l200ll : iLee et all l2009l : 
iTrenti et al.ll2010l ). Another possibility is that the galaxies in 
the "dark" halos are actively star-forming but are too dusty 
to satisfy LBG colour selections. In this case, the occupancy 
constrains the scatter in the dust column distribution. To 
see this, note that randomly rendering half of all actively 
star-forming galaxies too dusty to satisfy LBG selection cri- 
teria would shift the inferred HOD down by 50%. For suffi- 
ciently steep low- mass cutoffs, this would not change Afmin, 
hence the ratio Ng/Nh would also decrease by 50%. It is also 
possible that not all halos above Mmin host galaxies whose 
masses are comparable to those of LBGs (LBG "cousins"). 
In this case, {Ng)M constrains the CLF and a low observed 
occupancy has two possible interpretations: Either the CLF 
has a low normalization and many halos do not host galax- 
ies at all, or it has a steep faint end and "dark" halos host 
low-mass galaxies that are too faint to be observed. We will 
now demonstrate that our simulations naturally yield this 
last scenario, and that there is therefore no confiict between 
a low observed occupancy and a smoothly-rising SFH sce- 
nario. 

In order to illustrate this possibility, we shall consider 
our simulated LBG HOD sd, z — 5.5. To derive it, we identify 
dark matter halos using the friends-of-friends group finder 
FOIG with a linking length tied to the virial overdensity at 
z = 5.5. We then match the halo and galaxy catalogs and 
assign galaxies to parent halos. Finally, we count the num- 
ber of simulated galaxies brighter than absolute magnitude 
cuts of -18 and -19 that live within each halo. This simple 
approach does not isolate subhalos from parent halos, hence 
it only allows us to compute the sum of the two terms in 
Equation [S] However, this is sufficient for our present pur- 
poses; we leave a more detailed analysis of our simulated 
HOD for future work. 

In the top panel of Figure [13] we show the resulting pre- 
dicted HODs. Broadly, the simulated LBGs live in halos with 
masses in the range log{Mh/MQ) =10-12. Brighter galaxies 
live preferentially in more massive halos as observat ionally 
inferred (for example. ICiavalisco fc Dickinsonll2001^ . Both 
HODs rise steeply from a luminosity-dependent minimum 
halo mass Mmin, then turn over by log (M^/M©) = 11, 
then rise as a power law to higher masses. Mmin, which 
we define as the mass of the least massive halo that hosts 
at least one galaxy satisfying the luminosity threshold, is 
log(Mmin/M0) =10.0 and 10.5 for M1350 < -18 and -19, 
respectively, and the corresponding power law slopes are 
0.8 ±0.1 and 1.0 ± 0.2. 

The fact that our HOD lies well below unity at 
low masses has two important implications. First, our 
simulated HOD does not resemble conventional "central- 
satellite" models in which nearly all hal os more massive than 
Mmin host Js 1 galaxy (for example, iBerlind et al.l l2003l : 
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iKravtsov etail |2004| : IZheng et all 120051 ) because roughly 
half of the simulated LBGs satisfying current detection lim- 
its live in halos less massive than the mass scale at which 
the HOD crosses unity. Second, it implies that the mean 
occupancy (Equation |3]) is also less than unity. We confirm 
this in the bottom panel, which gives the cumulative occu- 
pancy at each halo mass for all halos between Mmin and that 
mass. Our predicted occupancies of 0.15-0.3 are consistent 
with the observation ally inferred range at 2 = 4-5 of 0.15- 
0.60 IjLee et al.ll2009f) . However, they do not correspond to a 
star formation duty cycle because our simulated duty cycle 
is essentially unity (Figures [5] and [TU]). Instead, they indi- 
cate that our simulated CLF does not assign an observable 
galaxy to each halo above Afmin. This possibility reconciles 
the smoothly-rising SFH scenario with the observed cluster- 
ing behavior of LBGs, but it raises the question of why our 
predicted HODs climb so much more slowly than would be 
expected in conventional central-satellite scenarios. 

The answer lies in the impact of outflows. DFO06 
showed that, in the absence of outflows, halos sd, z — 6 are 
expected to contain roughly their global fraction of baryons. 
These baryons invariably collapse into galaxies on a dy- 
namical timescale, leading t o the conventional HO D picture 
such as the SPH model of iBerlind et al.l (|2003l ). By con- 
trast, momentum-driven winds evacuate an increasing frac- 
tion of the baryons at lower halo masses because the mass- 
loading factor Tjw scales inversely with the velocity disper- 
sion. This mass dependence boosts the scatter in the rela- 
tionship between baryonic mass and halo mass because the 
instantaneous outflow strength depends on both the galaxy 
mass and the redshift, with the result that a galaxy's time- 
averaged outflow strength depends on its entire growth his- 
tory. Consequently, a significant fraction of isolated central 
galaxies that would be observable in the absence of outflows 
are suppressed below observational limits, rounding off the 
low-mass end of the HOD and suppressing the predicted 
occupancy. Broadly, these effects are expected to be signifi- 
cant if the outfiow mass- loading factor rjw > 1. For galaxies 
brighter than Muv = —18, our simulation assigns mass- 
loading factors of 3-7, hence this condition is satisfied for all 
observable galaxies at z — 5.5. Note that the possibility that 
feedback coul d dramatically alter the predicted HOD wa s 
also raised bv lBerlind et al.l(|2003h and lzheng et all l|2005h . 

While we have argued that scatter in halo formation his- 
tories dominates the halo occupancy at low halo masses, two 
other processes could scat ter galaxies below A/iaso — —18. 
First. lYoshida et al.l (120021 ) have shown that, in the absence 
of outflows, halos consisting of fewer than 75 dark matter 
particles show significant scatter in their gas fractions ow- 
ing to numerical effects. Such an effect would suppress the 
predicted occupancy at low halo masses. This effect should 
be negligible in our simulations because the smallest halo 
that hosts an observable galaxy contains roughly 1000 dark 
matter particles, indicating that our resolution cut is conser- 
vative. Second, intrinsic scatter in the relationship between 
dust column and halo mass has a similar effect because it se- 
lectively suppresses the luminosities of the central galaxies of 
some low-mass halos below detection thresholds. Our dust 
prescription does contain intrinsic scatter, and the HODs 
in Figure [13] include this effect. Neglecting dust extinction 
boosts the predicted occupancies for the fainter sample from 
0.3 up to 0.4. Hence the predominant factor suppressing the 
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Figure 14. The simulated CLF at 2 = 5.5 in three bins of halo 
mass; mean halo masses and occupancies within these bins are 
indicated. The CLF is well-approximated by a power law whose 
slope steepens and normalization decreases to lower masses. The 
mean occupancy increases with increasing halo mass as observed, 
but in our simulations this reflects the mass-dependence of the 
CLF rather than a star formation duty cycle. 



predicted occupancy of low-mass halos is outflows rather 
than numerical effects or dust extinction. 



7.2 Conditional Luminosity Function 

If "dark" halos do not host LBGs that are experiencing a 
temporarily quiescent phase, then what sort of galaxies do 
they host? We propose that halos that are more massive 
than Mjnin but that do not host an observable galaxy pref- 
erentially host galaxies that are simply faint. In Figure [HI 
we illustrate this possibility through the predicted CLF at 
z = 5.5 in three bins of halo mass. The halo mass bins are 
1.26 < Mh < 2.5, 2.5 < M^ < 5, and 5 < M^ in units 
of lO^^M©; the corresponding mean halo masses are indi- 
cated in the Figure. For each bin, we compute the predicted 
CLF by counting the number of galaxies as a function of 
absolute UV magnitude (including dust) whose host halo 
mass falls within that mass bin and then dividing by the 
number of simulated halos in that bin. Our CLF is well- 
approximated by a power law whose slope fiattens and nor- 
malization increases with increas i ng ha lo mass. Unlike the 
wind-free models of IZheng et al.l l|2005l ). it does not show 
a bump at high luminosities corresponding to the central 
galaxy contribution because the scatter in the relation be- 
tween halo mass and UV luminosity smooths this feature 
out. These CLFs have several interesting implications. First, 
their steep slopes imply that most halos above IO'^'^Mq host 
galaxies whose luminosities fall below current detection lim- 
its. Second, the presence of faint galaxies in massive halos 
(largely as satellites) explains the observation of enhanced 
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subst ructure around bright galaxies at « = 4 IjLee et al.1 
120061 . Figure 5). Finally, integrating these functions down 
to the current observational limit of M1350 <= —18.0 yields 
the mean observable occupancy in bins of halo mass, which 
we find grows from 0.05 to 1.56 as the halo mass increases 
by a factor of 10. This finding supports the assumption of 
HODs t hat predict a lo wer LBG occupancy at lower halo 
masses IjLee et al.ll2006l ). but within our simulations it re- 
flects the mass dependence of the CLF rather than a halo 
mass-dependent star formation duty cycle. 

In summary. Figures [13] and [TJ] show that our simu- 
lations are consistent with observed LBG occupancies of 
less than unity, but — given our near-unity star formation 
duty cycles — they predict that the "dark" halos preferen- 
tially host galaxies that are simply too faint to be observed 
rather than galaxies that are temporarily dusty or quies- 
cent. Of course, scatter in the dust column and SFR at a 
given halo mass must contribute to the observed occupancy 
at some level, but our simulations suggest that these effects 
are not dominant. It is possible that ot her feedback pro- 
ccsscs such as local ionizing backgrounds (|Cantalupdl2010l : 
[Gnedin 2010) or energy-driven outflows may be similarly 
successful in reproducing the observed occupancy, although 
momentum-driven outflows are attractive because they have 
already been shown to reproduce a wide variety of comple- 
mentary constraints (Section[2]). The important point is that 
there is no conflict between the observed low occupancies 
and the smoothly-rising SFH scenario. This interpretation 
makes several predictions, which we discuss in the next sec- 
tion. 



7.3 Predictions 

First, invoking strong outflows in order to explain the sup- 
pressed observed occupancy predicts that the minimum 
mass halo that can contribute to the observed sample Mmin 
is larger than would be expected given the galaxy's baryonic 
mass and the global baryon fraction. To illustrate this, we 
have determined the predicted relationship between M1350 
and baryonic mass at z — 5.5, multiplied by nji//f7(, to ob- 
tain the minimum mass halo that could host a galaxy with 
-^^i350=-18 and -19, and marked these minimum halo masses 
with arrows in the top panel Figure 1131 The actual mini- 
mum halo mass is roughly twice as large as would be ex- 
pected without strong fe edback (for example, iBerlind et al.l 
I2OO3I : IZheng et al.l [2005J) . Note that increasing our outflow 
strengths in order to improve the agreement with the ob- 
served colours and UV LF (Figure^ and [7]) would widen 
this gap. 

Second, mass-dependent outflows that "round off' the 
HOD at low masses tend to increase the minimum mass 
halo that hosts an average of one satellite galaxy, A/i , even 
more than they increase Mmin- Consequently, we find that 
the rat io Mi /Mmin i s 17- 18 for both simulated LBG sam- 
ples. IZheng et al.l l|2005h noted that simulations that do 
not treat outflows predict a ratio of 14 whereas a SAM that 
did predicted ratios of roughly 18. By contrast, low red- 
shift observations indicate a ratio of 18 (IZehavi et al.ll2005l : 
IZheng et al.ll2007l ). in excellent agreement with our outflow 
model. These considerations suggest that galaxy clustering 
measurements are sensitive to the effects of galactic outflows. 



although further inquiry is required in order to understand 
how constraining these measurements are. 

Third, if the galaxies living in "dark" halos are LBG 
"cousins" that happen to be quiescent or dusty, then 
the complementary galaxy fraction should be visible in 
optically-selected samples. Hence the smoothly-rising SFH 
scenario predicts that UV-select ed and optically-selected 
samples should largely overlap. iBrammer fc van DokkumI 
(,2007. ) used a Balmer break selection to identify samples 
of galaxies at 2; ~ 2.4 and 3.7 that directly test this idea. 
They found that, while some of their z ~ 2.4 sample possess 
red UV continua indicative of significant dust or a less active 
star formation phase, the majority of their z ~ 3.7 sample 
possess blue continua and would satisfy LBG colour cuts 
even though they were selected on the basis of their Balmer 
breaks. This finding argues against the existence of a signif- 
icant po pulation of quiescent or dusty B-dropout "cousins" 
(see also lOverzier et al.l [2009|) . Note that this conclusion is 
not in conflict with the observation that ^ 80% of massive 
(M* > 5 X 10^° Mq) galaxies at z ^ 3.5 may be too red for 
LBG selections (|Mancini et al.l 120091 ) because the majority 
of star-forming galaxies at z ^ 4 are not expected to be 
this massive-indeed, observations suggest that galaxies be- 
come increasin gly low-mass and dust-free at redshifts above 
3 (for example, Brammer fc van Dokkumll2007l : IVerma et al.l 



I2OO7I : iBouwens et al.ll2009al . l2010al ). Further exploration of 
optically-selected sam ples at z ^ 4 would be use f ul in c on- 
firming the results of iBrammer fc van DokkumI l|2007r ) as 
their samples were rather small (23 candidates at z ~ 4) 
and lacked spectroscopic redshifts. 

Fourth, the scatter in the SFR-M* relation should be 
small since short duty cycles would give rise to a significant 
number of galaxies that were transitioning between active 
and quiescent states. For example, assuming that galaxies 
transition between active and quiescent states on a halo 
dynamical time of « 100 Myr and that the ac tive star- 
formation period lasts 400 Myr (|Lee et al.ll2009l ). roughly 
half the number of actively star-forming LBGs should be ob- 
servable in a transition s t ate below the star-forming "main 
sequence". JLabbe et al.| (|2009j) measure a small scatter of 
0.25 dex at z = 7 and do not observe a transition popula- 
tion even though their observations are sensitive to galaxies 
well below the star-forming main sequence (for example, see 
the inset panel in their Figure 2). These observations are 
difficult to reconcile with short star formation duty cycles. 

Finally, we predict that the observed occupancy should 
not be lower for fainter samples than for brighter ones be- 
cause a larger fraction of halos host faint galaxies than bright 
galaxies (Figure [14]). This is visible in the bottom panel of 
Figure [13] where thin and thick curves indicate the cumula- 
tive occupancy without and with the contribution of satellite 
galaxies. If satellite galaxies are ignored then the occupancy, 
defined in this case as the fraction of halos above Afmin that 
host at least one observable galaxy, is roughly the same for 
fainter and brighter samples. If satellite galaxies are included 
then the occupancy, now defined as Ng/Nh{Mh > Mmin), 
increases for fainter samples owing to the contribution of 
satellites. Note that this prediction is in direct conflict with 
a scenario in which low-mass galaxies at high redshift pos- 
sess burs ty SFHs reminisc ent of local dwarf galaxies. Unfor- 
tunately, [Lee et al.l l|2009r ) assumed a single duty cycle and 
did not explore this question in detail. Intriguingly, how- 
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ever, their brightest subsample at 2 ~ 4 prefers duty cycles 
of 10-50% whereas their full sample — whose clustering prop- 
erties are presumably dominated by the more numerous faint 
galaxies — prefers longer duty cycles of 25-100% (their Fig- 
ure 9), in qualitative agreement with our prediction. More 
detailed exploration of the luminosity dependence of the 
occupancy should be possible with current and upcoming 
data sets and would constitute a useful test of our predicted 
SFHs. 



8 TRENDS AT LOWER REDSHIFTS 

Observations at z < 6 indicate that galaxies eventually de- 
part from the smoothly-rising, scale-invariant SFHs that are 
observed at earlier epochs. In this Section, we comment on 
these observations and discuss relevant physical processes 
that could emerge at 2 < 6. 

We begin by considering whether star formation is a 
steady-state process at all masses. The SFHs of actively 
star-forming galaxies more massive than 10^ Mq are likely 
to be predominantly smooth at all redshifts. Observations 
in support of this view include the tight observed scat- 
te r in the SFR-M« relation (< 0.3 dex) do wn to 2 = 
(JNoeske et al.ll2007al : [Brinchmann et al.ll2004 ) and the fact 
that visibly merging galaxies account fo r less than 30% of 
all star formation durin g 2 ~0. 24-0. 80 l| Jogee et al.l l2009l : 
see also lWolf et aLll2005h . 

At lower masses, evidence that dwarf galaxies experi- 
ence short-lived starbursts has historically promp ted many 
to fav or generally bursty SFHs for dwarfs (jSearle fc Sargenti 
1 19721 ). Starbursts could owe to interactions or to an ir- 
regular gas inflow rate. However, more recent studies em- 
ploying larger samples have thrown the significan ce of star- 
bursts in dwarfs' SFHs into doubt. For example, iLee et al.l 
l|2009r ) found that dwarfs whose emission lines possess high 
equivalent widths live preferentially in underdense regions. 
This is inconsistent with the view that interaction-induced 
starbursts dominate dwarfs' SFHs because interactions are 
more frequent in overdense regions. It suggests instead that 
star formation is suppressed in dwarfs that become satel- 
lite galaxies. This is expected if inflows drive star forma- 
tion because inflows funnel gas preferentially onto the cen- 
tral galaxy. Note that we expect inflows to be smooth at 
2 > 6 because most observed galaxies are central galaxies 
(compare thin and thick curves from the bottom panel of 
Figure [13). 

More recently, iLee et al.l l|2009l ) found that only 6% of 
a volume-complete sample of « 300 local dwarf galaxies are 
currently undergoing starbursts. Further, bursts accounted 
for at most one quarter of t he total current sta r formation in 
their sample. Together with ljogee et al.l (120091 ) . these results 
indicate that star formation in the central galaxies of ha- 
los below 10^^ Mq is a predominantly smooth process at all 
redshifts, which is a prediction of our model as well as oth - 
ers (for example, [Guo fc WhitdboOSi : ICattaneo et al.ll2010^ . 

Departures from rising SFHs are observed after 2 = 2 
in two senses. First, the emergen ce of a bimodality in 
galaxy colours as early as 2 = 2 (jGiallongo et al.l |2005| ) 
reflects the fact that many galaxies eventually cease form- 
ing stars. These " quenched" galaxies live preferent i ally i n 
overdense regions IjHogg et al.l |2004| : ICooper et al.l |2010|), 



which supports the view that quenching owes to the on- 
set of ho t-mode accretion i n halo s more massive than ~ 
10^^ Mq (|Keres et al.1 12005| . I2009D . po ssibly coupled with 
extra heating from a central source (|Croton et al.l I2OO6I : 
ICattaneo et al.ll2006l ). These effects should be weak in cur- 
rent 2 > 6 samples because their host halos are generally 
less massive than 10^^ Mq (Figure [13} • 

Second, observations indicate that the SFRs of galax- 
ies whose SFHs are n o t yet quenched decrease fo r 2 < 

1 (JNoeske et al.1 l2007al : JPerez-Gonzalez et al.1 I2OO8I '). This 
cannot be attributed to evolution in the interaction rates 
because obser vations indicate tha t the merger r ate evolves 
either mildly IjCassata et al.ll2005l ) or not at all IjLotz et al.l 
1200a ) out to 2 ~ 1. On the other hand, it occurs naturally 
if the gas accretion rate is assumed to trace the host halo's 
growth rate because halo growth rates fall off at late times 
owing to cosmological expansion. Models that incorporate 
this effect readily yield risi ng SFHs at early times and de- 
caying SFHs at lat e times (|Davell2008l : iBouche et"aLll2009l : 
iDutton et aLll2009^ . 

Finally, observations indicate that the slope of the 
SFR-Af, relation flattens f rom ~unity at 2 = 6 t o 0.7- 
0.9 at 2 < 2 (for example. iBrinchmann et al.ll2004h . This 
dep arture from scale-i nvariance is equivalent to downsiz- 
ing llCowie et al.||l996|). O ur simulations predict a slope of 
Ri 0.9 at 2 ^ 2 (Dave 2008), hence they qualitatively repro- 
duce the observed evolution. Two effects contribute. First, 
the fraction of galaxies living in massive (> 10^'^Mq) ha- 
los increases with time and mass. If the early stages of 
quenching (see above) involve suppressed inflow rates and 
a gradual exhaustion of available gas reservoirs, then the 
amount of suppression should increase with time and mass. 
This should mani fest both a s an overall flattening of the 
SFR-M* relation (lDavell2008r ). and as a dependence of the 
slope on environment. Second, feedback processes that de- 
lay gas accretion or star formation in low-mass halos flat- 
ten t he predicted slope (JNoeske et al.l l2007bl : JBouche et al.l 
[2OO9I). Momentum-driven outflows are a natural candidate 
because the mass-loading factors scale inversely with the ve- 
locity dispersion. This in turn causes gas in low-mass halos 
to spend most of its time suspended above the galaxy's ISM. 
Outflows are also appealing because they reconcile our sim- 
ulations wit h a wide rang e of co mplementary observations 
(Section [TJ. [Putton et al.l l|2009r ) argued that outflows can- 
not change the slope of the predicted SSFR trend, but their 
model did not treat the possibility of outflows that do not 
escape from the halo. 

While our simulations qualitatively reproduce the ob- 
served flattening trend, in detail they predict slightly less 
flattening than is observed. This could sugges t either that 
inflow s in massive halos are undersuppressed IjGabor et al.l 
[2OI0I), or that star formation in low-mass halos is delayed 
by other processes in addition to outflows. One possibil- 
ity is h eating by the local ionizing ba ckground within each 
galaxy l|Cantalupoll2010l : lGnedinll2010r ) . Our model does not 
treat the local background, but it would be interesting to 
explore in future work. 

In summary, star formation in halos below W^^Mq is 
likely to be a smooth process at all redshifts, but below 

2 «2-3 the SFRs begin to dechne and the SFR-M* trend 
flattens. The emergence of declining SFRs following 2 « 
2 owes to the decline in halo growth rates, which in turn 
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owes to cosmological expansion. Meanwhile, the departure 
from scale-invariance owes to at least two effects. In massive 
halos, accretion shocks and AGN feedback likely suppress 
inflows. In low-mass halos, outflows expel cold, star-forming 
gas while ionizing backgrounds heat it more rapidly than it 
can be replaced by fresh inflows. 



9 SUMMARY 

We have conducted a detailed comparison between the pre- 
dictions of state-of-the-art cosmological hydrodynamic simu- 
lations including momentum-driven outflows and recent ob- 
servations of reionization-epoch galaxies in order test two 
fundamental predictions regarding high-redshift SFHs; 

(i) SFHs at early times are smoothly-rising; and 
(ii) SFH shapes are scale-invariant. 

Our findings are as follows: 

(i) Our simulations reproduce the observed rest-frame 
UV luminosity function as well as its integral to current 
limits to within a factor of 2-3 during the epoch z =6-8. 
They also reproduce the observed stellar mass density to 
current limits. Either of these constraints could in princi- 
ple be matched by any SFH at a single epoch, but the only 
smooth SFH that can account for their evolution is a ris- 
ing one. This suggests that our prescriptions for star for- 
mation and feedback yield SFHs in reasonable agreement 
with current constraints. In detail, we find that increasing 
the assumed outflow mass-loading factors by «50% could 
bring the predicted and observed UV luminosity densities 
into better agreement without compromising the good agree- 
ment with the observed stellar mass density (to within the 
reported errors). Integrating below current detection limits, 
we estimate that current observations probe the brightest 
30% of the total star formation and stellar mass densities at 
z^6. 

(ii) We reproduce the observed blue UV continua at 2 = 7 
for all but the faintest galaxies and marginally reproduce the 
observed bluening to lower luminosities. The UV continua 
of the faintest galaxies may sugg est lower metallicities a nd a 
nonzero ionizing escape fraction (jBouwens et al.ll2010al ). We 
attribute the observed colour-magnitude trend to a modest 
level of dust redde ning that increases with luminosity as see n 
at lower redshifts (iMeurer et al.lll999l : [Shaplev et al.ll200l[ ). 
This suggests that exotic stellar populations are not required 
in order to explain observations of galaxies brighter than 
Hi60 = 27.5 at 2; = 7 (JFinkelstein et al.ll2009l : lBouwens et all 
l2010al ). 

(iii) Our simulations' relatively evolved stellar popula- 
tions predict red Hieo — [3.6] colours at z > 6, augmented 
by dust and optical emission lines. However, t he predicted 
colour s are ~ 0.5 mag bluer than reported bv lLabbe et al.l 
l|2009f ). It is not clear whether resolving this discrepancy re- 
quires significantly higher mass resolution or improved ob- 
servations. For the brightest galaxies, allowing the redshift 
and extinction to vary freely enables our models to repro- 
duce the observed -R'-[4.5] fluxes while underproducing the 
observed J125 and Hieo fluxes. This suggests that the red 
observed -ffieo ~ [3.6] colours may be as likely to owe to 
systematic error in the WFC3 as in the IRAC observations. 



(iv) We reproduce the near-unity slope of the observed 
SFR-M* relation at all redshifts. This agreement supports 
a scenario in which observed reionization-epoch galaxies be- 
gan forming stars at similar epochs and possess SFHs that 
vary only by a scale factor, which is a robust prediction 
of hydrodynamic simulations. We also reproduce the ob- 
served slow evolution of the SSFR, which is incompatible 
with constant or smoothl y declining SFHs while support- 
ing smoothly-rising SFHs (jstark et al][2009 '). The predicted 
scatter is ~ O.ldex, less than the (still remarkably tight) 
reported scatter of 0.25 dex. The larger observed scatter 
may imply a need for higher mass resolution in order to 
resolve the impact of minor mergers, but it may also be 
boosted by observational uncertainties. An offset of a factor 
of 2-3 between observations and predictions suggests that 
observationally-inferred stellar masses should be corrected 
for the likely impact of optical emission lines. 

(v) Our significantly improved treatment for the produc- 
tion and transport of metals predicts that observable galax- 
ies at z Js 6 possess metallicities in excess of O.IZq. This is 
several times larger than our previous prediction. Increas- 
ing the assumed outfiow mass-loading factors to improve 
the agreement between the predicted and observed UV LFs 
would lower these predictions by ~ 0.1 dex, but it would 
not change the fundamental result that galaxies self-enrich 
quite quickly. Hence we confirm that it should be possible 
to understand current samples without reference to Popu- 
lation III stars. Nevertheless, a residual fraction 10"'^ of all 
star formation a,t z = 6-8 may occur below the Population 
III metallicity threshold, which may eventually be visible as 
pair-production supernovae or GRBs. 

(vi) There is no conflict between observations suggesting 
that many massive halos do not contribute to the observed 
samples at z = 4-5 and the smoothly-rising SFH scenario. 
This is because clustering observations constrain the halo 
occupancy rather than the star formation duty cycle. The 
latter is only one of several possible interpretations of an 
overall occupancy that is less than unity, and our simula- 
tions disfavor it. Instead, we favor the interpretation that 
it reflects the scatter in the luminosity-halo mass relation. 
This is because strong outflows boost the scatter in the SFRs 
of low-mass halos, suppressing the UV luminosities of many 
isolated central galaxies that would be visible in the absence 
of outflows. Consequently, more than half of observed galax- 
ies live in halos whose mean (observable) occupancy is less 
than one and the overall fraction of halos above Mmin that 
host an observable galaxy is 0.2-0.3 as observed even though 
the predicted star formation duty cycle is unity. A halo oc- 
cupancy in this range should thus be regarded as indicative 
of a high star formation duty cycle rather than a low one. 
This interpretation makes a number of predictions: 

(a) The minimum mass for hosting an observed galaxy 
Mmin is at least twice as large as would be expected given 
galaxies' baryonic masses and the global baryon mass frac- 
tion; 

(b) The minimum mass for hosting an average of two 
galaxies is larger than Mmin by a factor of 17-18; 

(c) Optically-selected and UV-selected samples should 
largely overlap at 2 ^ 4; 

(d) The scatter in the SFR-M* relationship should be 
small, and few galaxies should be observed with sup- 
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pressed SSFRs indicative of a transition between active 
and quiescent states; 

(e) The overall occupancy for faint samples should be 
greater than or equal to the overall occupancy for brighter 
samples, in direct conflict with a picture in which lower- 
mass galaxies possess burstier SFHs in analogy with local 
dwarf galaxies. 

Each of these predictions is in tentative agreement with cur- 
rent observations, but we expect that more detailed analyses 
of current and upcoming samples will test them more thor- 
oughly. 

Further progress in understanding high-redshift SFHs 
will occur on both theoretical and observational fronts. 
On the theoretical side, it would be interesting to in- 
crease our mass resolution significantly while relaxing the 
assumption that wind particles temporarily decouple hydro- 
dynamically in order to resolve feedback on small scales. 
Simulations of dwarf galaxy formation at low redshift 
tend to find burstie r SFHs as the ma ss resolution is in- 
creased ffo r example. I Valcke et al.ll2008r) owing to stochastic 
processes (jGerola et al.l [l980|) , hence it is natural to sup- 
pose that similarly low-mass objects at z ^ 6 would also be 
bursty. On the other hand, dwarf galaxy SFRs are gener- 
ally < O.5M0yr~^ whereas even the faintest z = 7 galaxies 
indicate SFRs of > 2M0yr~^; whether bursty SFRs are ex- 
pected in the presence of efficient gas inflows is unclear. If 
so, this may improve the agreement with the observed scat- 
ter in the SF R-M« rela tion as well as the colours of low- 
mass galaxies (jLabbe et al. 2009). Unfortunately, increasing 
our mass resolution is not trivial because our present sim- 
ulations already probe current computational limitations. 
Further progress will require resimulations o f a statistical 
samp le of high-redshift galaxies (for example, iBrooks et al.l 
[20071). By contrast, incorporating stronger outflows in or- 
der to improve the agreement with the observed UV lumi- 
nosity function would be trivial (aside from the computa- 
tional expense), although it remains to be seen whether this 
would also preserve the agreement with the observe d density 
of met als in the intergalactic medium (jOppenheimer et al.l 
l2009bl ). 

The impact of an inhomogeneous ionizing background 
warrants closer scrutiny. Our curre nt work does not indicate 
that our assumption of a uniform iHaardt fc Madaul (|200l|) 
ionizing background produces large errors in our predictions 
(DFO06). Nonetheless, it would be useful to study the rela- 
tive roles of photoionization heating versus outflows in sup- 
pressing star formation in faint galaxies using simulations 
that treat the nascent ionizing background self-consistently. 
Preliminary work suggests that interesting nonlinear inter- 
actions exist between outflows and the ionizing background ; 
these warrant more detailed study (|Pawlik fc Schavell2009r ). 

A closer inquiry into the nature of our predicted HOD 
and its ability to reproduce clustering constraints at z = 4-5 
is in order. We have not delved into a full clustering analysis 
because our goal has been to demonstrate that clustering 
properties in the presence of mass-dependent outflows are 
qualitatively different than the conventional central-satellite 
HOD picture and have the potential to reconcile the ob- 
served duty cycle with smoothly-rising SFHs. To this end, we 
have already demonstrated encouraging agreement between 
our predictions and the observed occupancy and Mi/Mmin. 



A more detailed study of the CLF and its associated HOD 
in simulations with and without outflows would bring the 
impact of outflows on the HOD into sharper relief while 
providing a new class of template HODs for interpreting ob- 
servations. Such a study would involve using our predicted 
HOD to populate halos from a larger-volume N-body simu- 
lation in order to provide good statistics, a task that is well 
beyond the scope of our current work. 

On the observational side, further progress will follow 
from a more systematic reconstruction of the high-redshift 
SFR-M* correlation and its evolution, which are the most 
important observational probes of the "typical" SFH in an 
epoch where the SFHs of individual galaxies remain uncon- 
strained. This involves work on three fronts: First, deeper 
rest-frame optical constraints on larger samples of galax- 
ies are required in order to leverage existing UV-selected 
samples. Second, progress requires improved SED model- 
ing that invokes physically-motivated SFHs (rather than, 
for example, instantaneous bursts) and accounts for nebular 
emission in order to reduce systematic errors. Third, uni- 
form SED-fitting studies incorporating common modeling 
assumptions (for example, Schaerer & do Barros 2010) will 
reduce artificial scatter owing to differing modeling assump- 
tions. Because of the importance of the SFR-M* plot, we 
heartily encourage the publication of inferred physical prop- 
erties order to facilitate further theoretical comparisons. 

In order to gain a better understanding of the scat- 
ter in the SFR-M* relation as well as the implications 
of the low observed halo occupancy, we encourage fur- 
ther inquiry into the overlap between UV-selected and 
optica lly-selected samples at z ^ 4, following the exam- 
ple of iBrammer fc van DokkumI (|2007|) . These studies are 
a strong test of our prediction that SFHs are predominantly 
smooth rather than bursty. 

With the growing body of observational evidence that 
high-redshift SFHs are smooth rather than bursty and the 
implication that the progenitors of F-dropouts, for exam- 
ple, are visible as i-dropouts, it is important that SED- 
fitting studies account for galaxies' hypothetical progeni- 
tors at earlier epochs in addition to their observed fluxes. 
For example, exponentially decaying SFHs with a decay 
time of 300 Myr applied to V-dropouts predict bright pro- 
genitor i-dropouts whose colours/SSFRs are observationally 
ruled out (unless they are very dusty, which again reinforces 
the need for study of optically-selected samples). For the 
same reason, it is now becoming possible to constrain high- 
redshift galaxy SFHs using stati stical methods applied at 
multiple epochs (see, f or example, Conrov fc Wechslerll2009l : 
iMaraston et all |201G| : iPapovich et al.l I2OI0I '). Such studies 
provide complementary tests of our predicted SFHs. 

Finally, we have outlined a number of predictions re- 
garding the clustering properties of faint (M1350 > —20) 
high-redshift galaxies that may be amenable to study us- 
ing current samples (for example, from GOODS) and that 
should be considered in the design of upcoming Hubble and 
Spitzer programs to study the distant universe. 
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